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1. WPROWADZENIE 
1.1. Cel i plan pracy 

Opis układów niejednorodnych stanowi od lat jeden z głów- 
nych nurtów fizyki statystycznej płynów. Jest on kluczem: do 
zrozumienia zjawisk powierzchniowych i przejść fazowych zarówno 
w płynach prostych jak i układach ciekłokrystalicznych, 

W pracy zostanie sformułowana mikroskopowa teoria SDA 
/smoothed density approximation = przybliżenie wygładzonej 
gęstości/ dla opisu niejednorodnych i anizotropowych układów 
złożonych z twardych cząstek. W ramach tej teorii zbadamy 
przejścia fazowe w układzie twardych sferocylindrów oraz twar- 
dych elipsoid i wyjaśnimy mikroskopowe mechanizmy, które prowa- 
dzą do pojawienia się różnych faz ciekłokrystalicznych. W nisko- 
gęstościowej granicy Onsagera zbadamy jak i dlaczego cząsteczki 
ciekłego kryształu porządkują się przy powierzchni oraz obli- 
czymy napięcia powierzchniowe fazy nematycznej współistniejącej 
z faza izotropową oraz fazy nematycznej w kontakcie z twardą 
ścianą. Zbadamy także profile gęstości i parametrów uporządkowa- 
nia przy ścianie, zwracając uwagę na efekty nielokalne. 

Na zakończenie przedstawimy metodę umożliwiającą praktyczne 
stosowanie pełnej teorii SDA w badaniu zjawisk powierzchniowych 
w układach gęstych. 

Układ pracy jest następujący: 

W rozdziale pierwszym omówimy klasyfikację faz ciekłokrys- 
talicznych z podaniem ich symetrii. Rozdział drugi poświądóny 
jest netodon funkcjonałów gęstości w układach anizotropowych 
i niejednorodnych. VWyprowadzimy w nim formalne wyrażenia na 
wielki potencjał termodynamiczny jako funkcjonał jednocząstecz- 


kowej zredukowanej funkcji rozkładu. 
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Pokażemy, że warunek na minimum tego funkcjonału odpowiada 
warunkowi równowagi w układzie, oraz podamy warunki stabilności 
dla układu niejednorodnego i anizotropowego. Rozdział trzeci 
zawiera elementy fizyki statystycznej powierzchni, a w szcze- 
gólności wyprowadzenia ścisłych reguł sum dla płynu anizotro- 
powego w kontakcie z twardą ścianą. Analiza modelu Onsagera i 
teorii niejednorodnych płynów prostych doprowadzi w rozdziale 
czwartym do sformułowania teorii SDA dla niejednorodnych i ani- 
zotropowych układów złożonych z twardych cząstek. W rozdziale 
piątym zastosujemy SDA w granicy małych gęstości do badania 
zjawisk powierzchniowych w ciekłych kryształach. Uzyskane wyniki 
porównamy z wynikami doświadczeń, W rozdziale szóstym zbadamy 
przejścia fazowe w gęstym układzie twardych sferocylindrów i 
elipsoid obrotowych. Wyniki obliczeń porównamy z symulacjami 
komputerowymi. Dyskusję wyników pracy przedstawimy w rozdziale 


siódmym. 


1.2. Charakterystyka faz ciekłokrystalicznych 

Ciekłe kryształy są cieczami anizotropowymi. Składają się 
z wydłużonych /lub płaskich/ cząsteczek związków organicznych, 
które mogą się porządkować wzdłuż wyróżnionego kierunku w przes- 
trzeni. Ten typ uporządkowania charakterystyczny dla wszystkich 
ciekłych kryształów jest podstawą ich anizotropii. Wedle 
klasyfikacji wprowadzonej przez Friedla /1/ wyróżniamy trzy 
podstawowe typy faz ciekłokrystalicznych. Sa to nematyki, choles- 


a 


teryki i smektyki. 
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Nematyki /N/ 
W fazie nematycznej obserwujemy długozasięgowy porzedek orien- 


tacyjny /rys.l.l.,1.2./ przy jednoczesnym braku długozasięgowego 


porządku translacyjnego. 


RNA 


rys.l.l. faza nematyczna rys.1.2. faza nematyczna 
- cząstki wydłużone - Cząstki płaskie 


Cząsteczki układają się wzdłuż kierunku Å zwanego "direktoren". 
Oprócz osi symetrii rotacyjnej określonej przez ñ, układ posiada 
również płaszczyzny symetrii równoległe i prostopadłe do ñ. 
Kierunki Ń i -ñ są równoważne. Stopień uporządkowania w nematyku 


opisujemy parametrem porzedku /2/,Q, 
Q= < P(A6)) «JE (A-6)-F(A-O)dO, (MM) 


gdzie Po jest drugim wielomianem Legendre'a, () wersorem osi 
cząstki a f orientacyjfig funkcją rozkładu spełniającą warunek 


normalizacji 


OOE NE (4.2) 


Dla fazy izotropowej /I/ f=1/4TT i Q=0. 
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SE 


W ogólności aby opisać uporządkowanie orientacyjne niejedno- 
rodnej fazy nematycznej /np. przy twardej ścianie/ wprowadza 


się tensor parametru porządku /2/, 


Qj= ZŚ3wiuj - 84) (13) 


który w lokalnym układzie współrzędnych / All 2/ zapisujemy 


v postaci 


Qy=| 0  -z(8-P) O (1.4) 
0 o Q 


gdzie parametr porządku P, zwany dwuosiowością opisuje upo- 
rządkowanie orientacyjne cząsteczek w płaszczyźnie prostopadłej 
do A. 

Cholesteryki /Ch/ 

Lokalnie faza cholesteryczna przypomina fazę nematyczne, jednak 
kierunek Å zmienia się periodycznie w przestrzeni rys. /1.3/ 


według równania 
Mx = (qZ + 9) 
Ty = sin (qoz + 9) (4.5) 
ng =0 


Oś z jest tutaj osią skrętności /helical axis/. 
Okres przestrzenny L= T49 jest rzędu długości fali świetlnej 
X ~ 4000Å i jest około stu razy większy niż typowa długość 


cząsteczki, 
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rys.1.3 faza cholesteryczna 
„Molekuły tworzące fazę cholesteryczną nie mają symetrii 
zwierciadlanej. Cholesteryki sę interesującym. obiektem badań 
z punktu widzenia teorii defektów /3-6/. W fazie cholesterycznej 
obserwuje się spontaniczne uporządkowanie linii dysklinacji 
w sieć kubiczną prowadzące do powstanie faz błękitnych /6-8/. 
Fazy błękitne sę rzadkim przypadkiem stabilnego uporządkowania 
defektów. 
Smektyki /Sm/ 
Fazy smektyczne są strukturami warstwowymi i tym różnią się 


ód poprzednio omówionej fazy nematycznej /rys.l.4/. 


rys.1l.4 Smektyk. A 
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Grubość warstw, d, jest rzędu długości cząsteczki, W zależności 
od sposobu uporządkowania cząsteczek w warstwie wyróżniamy 

trzy główne typy smektyków: A,B,C, W smektyku A /SmA/ wersor ń 
jest prostopadły do warstw. Środki mas cząstek umieszczonych 

w warstwie nie wykazują długózasięgowego porządku translacyjnego. 
Warstwy mają strukturę dwuwymiarowego płynu i mogą się przesuwać 
względem siebie. W smektyku C /SmC/ rys.1.5/ wersor fi nie jest 
prostopadły do warstw i dlatego SmC ma symetry dwuosiową 


w przeciwieństwie do SmA/. 


rys.l.5 smektyk C 


W smektyku B /SmB/ w odróżnieniu od SmA i SmC środki mas cząs- 
tek w warstwie tworze, strukturę heksagonalną lub quasiheksago- 
nalną. W układach złożonych z cząstek płaskich zamiast fazy 
smektycznej obserwuje się fazy kolumnowe /9/. Cząsteczki w tej 
fazie układają się w równoległe kolumny, które tworzą dwuwymia— 
rową sieć przestrzenną. Wzdłuż kolumny środki mas cząsteczek 

są rozłożone w sposób przypadkowy. Do tej pory wykryto fazę 
kolumnową h, która posiada uporządkowanie heksagonalne oraz 

fazę r o porządku prostokątnym /10/. Dokładne omówienie symetrii 


różnych faz ciekłokrystalicznych można znaleźć w pracach /2,11/. 
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2. TEORIA STATYSTYCZNA FAZ CIEKŁOKRYSTALICZNYCH, METODA 
FUNKCIONAŁÓW. GĘSTOŚCI 


Teoria funkcjonałów gęstości /12-16/ /DFT - density functio- 
nal theory/ w fizyce statystycznej płynów prostych, pomyślana 
została jako metoda badania układów niejednorodnych, Wielkością, 
która określa stopień niejednorodności w układzie jest lokalna 
gęstość P .„ W ramach DFT /16/ wprowadza się dwa funkcjonały, 
Ay iF , lokalnej gęstości. Równanie na równowagowy profil 
gęstości gn go) otrzymujemy minimalizując (ly jeśli poten- 
cjał chemiczny jest ustalony lub minimalizując J jeśli ustalona 
jest liczba cząstek .w układzie ; (Ly [ fol] jest równy wielkiemu 
potencjałowi termodynamicznemu, a Fle] - energii swobodnej. 
DFT sformułowano również dla układów składających się z cząstek 
anizotropowych /17-24/ i temu zagadnieniu poświęcony jest roz- 
dział. 

W ogólności cząstka anizotropowa, nieliniowa, bez struktury 
wewnętrznej posiada 6 stopni swobody: trzy translacyjne i trzy 
rotacyjne. Jej położenie w przestrzeni określa wektor środka 
masy r i trży kąty Eulera 8,9,X - Ponieważ konkretne obliczenia 
przeprowadzamy dla układów cząstek o symetrii cylindrycznej /5 
stopni swobody, cząstka liniowa/, dla ustalenia uwagi ograniczymy 


formalne rozważania do tego przypadku. 


2.1. Podstawowe definicje. Gęstość prawdopodobieństwa i zreduko- 
wane funkcje rozkładu. 


Rozważamy klasyczny układ w równowadze, składający się z 
identycznych wydłużonych cząstek liniowych o masie m, zamknię- 
tych w objętości V i umieszczonych w polu zewnętrznym Vext* 
Mikroskopowy stan układu opisywany jest przez punkt w przestrze- 


ni fazowej r = (Nixe sexy) gdzie x= (£; 0, „p aBa Pa: Pe: ) 
określa położenieyjr,, pęd sprzężony kanonicznie z położeniem, 
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Pi orientację Oi: Pi oraz pedy sprzezone kanonicznie z NAA 
i-tej częstki. Zakładamy, że hamiltonian układu Hy (I) jest 
sumą trzech składników: energii kinetycznej, T, energii oddzia- 


ływania,U, potencjału zewnętrznego, Vis 


Hu=T+U+V ) (2.4) 
N 2 2 

LÅ + A+), (2.2) 
du $(5-4,0:,9:,6; gi) ) (2.3) 

Vi å Vere (TaSi), (2.4) 


gdzie Polti- Gi, 9:94) „jest dwuciałowym potencjałem oddzia- 
ływania, a I-momentem bezwładności cząstki względem osi przecho- 
dzącej przez środek masy cząstki i prostopadłej do jej osi syne- 
trii. Gęstość prawdopodobieństwa znalezienia układu w stanie mikro- 


skopowym P jest dane przez gęstość wiélkiego rozkładu kanoniczne» `` 


R (r)< SPGBH BAN , (2.5) 


ud 


go: 


scharakteryzowanego przez potencjał chemiczny Jų 1 temperaturę 
T =(kg Br gdzie kp jest stałą Boltzmana. Z} jest wielką sumą 


statystyczną określoną przez warunek normalizacji P,(T) : 
ferger) = ie fox... dry BUP) 4 (2.6) 


'dzie dx; =dr, dp; dQ dg; dp; dPy; a h jest stałą Plancka. 
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Wielkość P,(P) w pełni charakteryzuje układ; wartość średnia 


dowolnej wielkości mikroskopowej A (r) wyraża się następującym 


wzorem: 


<A» = fdr A(R) (23) 


Rozważmy A(P) w postaci sumy 4 
A(T)= ao +$ a406)+) On (XX). . . (2.8) 
L> tej=1 - 
Można sprawdzić ,że 


£ Ah på a |dX....dXk OG Ke) Fa (Xe), (2.9) 
gdzie 


Film M)=> NR EZ ~- -d Xy Rr) (240) 


dla k=0,1.... Niech we wzorze (2.9) Qr (--..Xk)E Qk (094... JG PK) n 
Wprowadźmy kat bryłowy w=(6,9) ) dw= ginOdedy z 


Wtedy ` 
- (2.44) 
SA», = a å far, day. .. digd QK (0 BOKO k (04 --- Ag); 
z 
gdzie k-cząstkowe zredukowane funkcje rozkładu Pk otrzymujemy 
przez całkowanie Fr PO zmiennych pędowych. Funkcje Ok posiada- 
rostą interprétac fi ną e.g. 1 
ja P a pretację zyczną e.g 24 (4 Gy) dydw 
jest prawdopodobieństwem znalezienia Śródka masy cząstki w ele- 
mencie objętości dra wokół r,, a jej osi symetrii w elemencie 
kata bryłowego duy wokół (4. Funkcję Pa(% 04) zapisujemy 
w postaci iloczynu lokalnej gęstości (r) i orientacyjnej 
funkcji rozkładu Z) spełniającej warunek normalizacji: 


faofa 
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W fazie izotropowej qlrw)= PT „ aw fazie nematycznej jedno- 
rodnej futro) = pf). W fazie smektycznej A (CAZ) jest funkcją 


periodyczną o okresie d: 
P1(1,W) "PA (qi+md,w), (2.42) 


gdzie m jest dowolny liczbą całkowitą, a fi "direktorem". 


Zdefiniujmy funkcje korelacji ha (170 0) następującym wzorem; 


h - man) _ 2.4 
NE aa | AE 


Wielkość ta opisuje odchylenie położenia par cząstek od roz- 
kładu przypadkowego i znika w granicy rE ESR „ W przybliże- 


niu niskogęstościowym 
h, =ePó_j=f,, (244) 


Funkcję f2(5-5,004) nazywany funkcją Mayera. Funkcja h, spełnia 
następującą regułę sum dla ściśliwości izotermicznej, Xr » /25/: 


-À x -A 99 = 
P (9%r)-8 (35)r (2.45) 
=Å+ N | dydey dv, de, hon a) ly) Qba, ), 


gdzie p jest ciśnieniem w układzie. 


2.2. Zasady wariacyjne dla funkcjonała wielkiego potencjału 
fermodynamicznago " || 
Rozważmy następujący funkcjonał: 

-A 
Q.[P] =<Hy=uN +B"lnP), (2.46) 


gdzie P jest dowolną gęstością prawdopodobieństwa spełniającą 
warunek normalizacji (2.6). Korzystając z definicji (2.5) otrzy- 


mujemy związek 
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QEPJE] +p lne} (2.47 


Ponieważ ng 0 dla PHPO więc QLP] przyjmuje minimalną 
o 


wartość dla PzP, równą wielkiemu potencjałowi termodynamicz- 


nemu, KOA 5 
Q=QIR]=-flnE, (248) 


Gęstość prawdopodobieństwa Pó (2.1) - (2.6) jest funkcjonałem 
Vexe" Łatwo sprawdzić, że zmiana Vext nawet o stałą prowadzi 
do zmiany P,. P, możemy również określić jako funkcjonał 
POE Uri) ponieważ g! jest jednoznacznie wyznaczona przez 
Vęxt /16,26/. Korzystając z tego faktu skontruujmy dwa funk- 
cjonały F " pe jednocząsteczkowej zredukowanej 
funkcji rozkładu (50) : 


FIgtma]=<T+U +88), (219) 


0,1(55)]= Fig] +|drdoglrojig (nu) „ujdrdogiru) (2.20) 


dla ustalonego V e Tak samo jak funkcjonat (2[P] przyjmuje 


ext 
wartość minimalng dla PsP, tak też Qv ma minimum dla 


fro) = QA( W) : 
SQv[(r0) 


jp 
„24 
Ó CO) Piro) Pirow) 2 ) 


Oy LE o= (2.22) 


Z definicji (2.19) i równania (2.21) wynika związek 


oraz 
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(ERESI pe Ven (moko . (223) 
pre) PT) 


Energia swobodna Helmholtza wyraża się wzorem 


F= fardo gu) Varm) TTP] . (2.24) 


Potencjał V występujący w równaniach (2.20) - (2.24) odgrywać 


ext 
może wiele ról. Właściwie dobrany symuluje ścianki naczynia; 
służy do opisu oddziaływania układu z dowolnymi powierzchniami 
ciała stałego m.in. twardą ścianą; lokalizuje powierzchnię roz- 
działu faz; przy badaniu przejść fazowych jednoznacznie określa 
rozwiązanie nieliniowych równań (2.21),£.23) .„ Stosowanie tych 
równań wymaga znajomości funkcjonału F . Zbadajmy bliżej 
jego strukturę, 


Dla U=0 rozważany układ redukuje się do idealnego gazu cząstek 


anizotropowych. Korzystając z (2.92.11) dostajemy Fo): 
PT) ap -B Vexe (r,0)) , (2 25) 
2 2 
N=(KE] (zę) (2.26) 


21m 211 


i zgodnie z definicja (2.19) otrzymujemy funkcjonał energii 
swobodnej Ha dla gazu idealnego 


Falstro]=p"fdrdo Aro an pru) - dj. (2.20 


Zdefiniujmy funkcjonał energii swobodnej związanej z oddziały- 


waniami w następujący sposób: 


F” [gro FIotzo)]-F tro],  ( 2.28) 


a także wprowadźmy funkcję 
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Przepiszmy równania (2.2]42.22) korzystając z tych wielkości, 


QF” Ana] +8 fardotw (u (wzięto) 4),(2.30) 


$(qo)=AexplBu BVy G w) + Gl m W; Loo] 3j s (2. 34) 


Funkcja c, ma interpretację efektywnego potencjału w polu 
którego znajduje się częstka w punkcie © o orientacji W. 
Potencjał ten jest wypadkową oddziaływania danej cząstki z po- 
zostałymi cząstkami układu. Z równania (2.31) wynika, że nawet 
jeśli Vęxt=0 to stan równowagi układu może odpowiadać konfigura- 
cji niejednorodnej dzięki efektywnemu potencjałowi ZE Jednak 
dla Vyxp=0 równanie to jest także spełnipne dla 9% =const. 
Aby znaleźć właściwe rozwiązanie (2.31) musimy albo wprowadzić 
pole łamiące symetrię /symmetry breaking field/ albo założyć 

z góry jakąś postać 91 „zależną od parametrów, o spodziewanej 
symetrii. W praktyce stosuje się tę drugą metodę, przy czym 

F lub Ny minimalizuje sie' względem tych parametrów. 


1.3 Funkeja korelacji Ornsteina-Zernike. Warunki stabilności 
az ciekłokrystalicznyc 


c4 jest pierwszą funkcją w hierarchii bezpośrednich funkcji 
ex Qx 
korelacji, Ch» generowanych przez 4. Różniczkując F n=krot- 


nie otrzymujemy 


Cn (10 „On; Leiro) "a al å nål 0) 2. (232) 
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Funkcja Sp (Gad ;L9YGCw)]) jest bezpośrednią funkcją korela- 
cji Ornstéina-Zernike /0-Z/. Różniczkując (2.31) względem 
g rw) dostajemy związek 


E-a) lor) _g Sulm, 
CC jay | eden), (233 


gdzie aa -Vxt* Z definicji (2.94-2.11) łatwo pokazać, że 


E Ua 
era Qala 410) + 993 605-%)6(0-0w)- 
= Cian) F lga) (2:34) 


co wraz z równaniem (2+33) oraz tożsamością 


Jada o ue =S(-BSlara) (2.35 
sukket dy GSA GaGa id 
ha Caen) = Cz (ey nød) + 

t Jorde, h (hett) (go) (236) 


wiążącego funkcję korelacji h (2.13) z c dla anizotropowego i 
niejednorodnego płynu. Korzystając z (2.36) przekształómy 
(2.15) dla fazy izotropowej tak aby ha zastąpić przez cz. 
W wyniku otrzymujemy związek pomiędzy Å ai a ściśliwością 


w fazie izotropowej /20/ 
Ble) BRB) =4— ojc (ra, a) dra døde (23) 
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Dodatniość prawej strony (2.37) określa stabilność mechaniczną 
fazy izotropowej. Poniżej podamy ogólny związek c Z warunkami 
stabilności układu ciekłokrystalicznego. 

w 8 2.2. pokazaliśmy, że dla V,,,=0 rozwigzanie równania (2.21) 
nie jest jednoznaczne. W szczególności aby danemu rozwiązaniu 
równania (2.21) odpowiadało minimum {Ly „ macierz drugich 
pochodnych Qy musi być dodatnio określona. Oznacza to, że 


nierówność 


fan duy dz da, Spion) Spk. 0») EE ez > O (2 38) 


musi być spełniona dla dowolnych zaburzeń Eoln) . Po przekształ- 
ceniu (2.38) /patrz (2.20), (2.27), i (2.32)/ otrzymujeny 


- związek 
Ja day dz do, S9C5w)8PL5)E Spion) 3 


- (7%, [Crop] >O. (2.39) 


Jest to warunek stabilności dla układu ciekłokrystalicznego. 
Jeżeli dla określonego potencjału chemicznego U= Up nierówność 
(2.39) przechodzi w równość, oznacza to, że faza ciekłokrysta- 
liczna opisywana przez omu) staje się niestabilna względem 
zaburzeń Soro). Wynika stąd, że dla Sil, musi nastąpić 
przejście fazowe od fazy o symetrii ru) do fazy o symetrii 
59 (Gu) „ Jeżeli przejście fazowe jest ciągłe to określa je 
warunek [L=/ł, . Analizę rozwiązania Piw) w pobliżu jj, nazywany 
analizę bifurkacji. Warunki stabilności dla układu ciekło- 


krystalicznego zostały sformułowane po raz pierwszy w pracach 
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I. Steckiego i A. Kloczkowskiego /18,19,20,27/. 


2.4. Gęstość funkcjonału energii swobodnej 
: ax > 
Załóżmy, że znamy wartość F” dla Prw). Całkując 
funkcjonalnie równość (2.29) wzdłuż pewnej drogi od Of W) 
x 
do Póriw) otrzynujenyJ Č pr] - Ponieważ F° jest jedno- 
znacznym funkcjonałem (iu) . to jego wartość w Pio) nie zależy 


od drogi całkowania. Wybierzmy więc drogę liniową 
Glod) = SY) + (Piriw)-PYGGu)), Ode. (2.40) 


W wyniku otrzymujemy 


Fegco]=FTĘ"xo]- faa Jagati- 
= $A) (00, [$05034)]). (2.44) 


Ta sama procedura prowadzi do związku /patrz (2.29) i (2.32) 

(Gz 9ir0d)])=cClGwiP0)]) + 

| add |dgdaj | glo, )- PY(G w, fx (2.42) 
KC (608426 51 90G0ydd')]), 


Stąd 
IF [prrd)]= F*TP(w)]-B" Jardo] PW) a PAlW)Ę x 
*Q(GETCWJ) —B Jdgdw dzda f oa-r 
P | 
følgen) på Bi) j dsa) del G (10%, Gd [lgay a'l), 
ò ò 


(2.43) 
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Spróbujmy teraz zdefiniować gęstość funkcjonału energii swo- 


bodnej, APT; Ego). równaniem 


FT gt] =foarers)avtr; eo). (2.44) 


Z równania (2.43) wynika, że funkcjonał AV jest nielokalny 

w skali zasięgu Ca» Załóżmy, że c ma zasięg potencjału 
oddziaływania Og . wtedy AP(x;[glxi0)]) zależy od (rw) dla 
|r-r|<Q. gdzie 2 jest zasięgiem $, e Oczywiście w fazie upo- 
rządkowanej możemy się spodziewać korelacji rozciągających się 
na obszarze rzędu rozmiarów układu, Fakt ten pominiemy przy 
konstruowaniu, w rozdziale 4., przybliżonego funkcjonału 
ILSA). 

W sposób równoważny formułuje się metodę funkcjonałów gęstości 
w zespole kanonicznym ,który w sensie średnich jest równoważny 
wielkiemu zespołowi kanonicznemu. W zespole kanonicznym F 
jest funkcjonałem, który minimalizujemy względem gixa) przy 
ustalonej średniej gęstości cząstek w układzie. Otrzymujemy 
wtedy równanie (2.23), w którym /Ł jest funkcją gęstości 
MG) . Warunek stabilności (2.39) również pozostaje bez 
zmian /18/. 
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3. ELEMENTY FIZYKI STATYSTYCZNEJ POWIERZCHNI W UKŁADACH 
ZŁOŻONYCH Z ANIZOTROPOWYCH CZĄSTEK 


3.1. Tensor napięć i równanie hydrostatyki 
Przypomnijmy wpierw znane definicje i zależności /28,29/. 
Mikroskopowy operator prądu cząstek, J(rit|M) , dany jest 


Jirel) => p Ś(0-1i(b)), (3.4) 


gdzie Pi jest pedem i-tej czastki a ri jej położenien. 
nn N 
Zniana J w czasie, uśredniona po przestrzeni fazowej [" z rozkła- 


dem P (2.5) Ke równa zeru, 


( 2x I =0 i (32) 


ponieważ w stanie równowagi nie może być makroskopowych prze- 


pływów. Z drugiej strony możemy wykonać formalny rachunek 


L om 
Ay (2 3 Sna) - 
N 
- 2 Er Gr, 


gdzie operator Ve jest gradienten w kierunku r, a 


O SE ro) = = > Ve (42, Å, (1-5, 0407) + 
+ Vee (0x) ) ))6 o Y OD 


(3.3) 


(34) 


Równanie (3.2) przyjmuje postać: 
Z. AA N 
c Å Gloop) sen)» 3 BL 


PV SC GP — Gr fr: ex uwi or) dlr-r)»=0 


Od 


(35) 


- 23 =~ 


Pierwszy człon równania (3.5) zapisujemy jako dywergencje 
tensora drugiego rzędu Po (r). /29,30/ 


VE =, 42 PA KAZ WW DE 
fat ŚCE-4)4, ) 


å ; 
który nazywamy konfiguracyjną częścią tensora napięć. C14 


(3.6) 


jest dowolnym konturem łączącym (i Z Eje 


Drugi człon w równaniu (3.5) jest dywergencją kinetycznej 


części tensora ciśnień /napięć/ Pulz)» zdefiniowanego wzorem 
is LODE x 
Pdr) 2 BER sitt, (33) 


gdzie © oznacza iloczyn tensorowy. Suma P, i P, definiuje 
tensor napięć p. Wstawiając (3.6)43.7) do równania (3.5) otrzy- 


mujemy równanie hydrostatyki /31/. 


Ve Plz) =- Jaro PA, 3) Ver Vare (710) (3.8) 


Tensor napięć E(r) nie jest wielkością jednoznacznie określoną 
ponieważ zależy od wyboru konturu całkowania Cij* Równanie 
hydrostatyki wyznacza jednoznacznie tylko jego dywergencję. 

W odróżnieniu od tensora napięć w płynach prostych, plr) w pły- 
nach anizotropowych może być niesymetryczny. 

Zarówno Ec jak i Ek wyrażają się przez zredukowane funkcje 


rozkładu /30./ 
R. 
R=- 7 [ana |do faw, Vi, Pa C0) x 
X 
Jet alr y T-L + fia 04100, ) | 
2 


P) = kT) T , (3.40) 
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SZĄ. 


gdzie I jest tensorem jednostkowym. Dla układu jednorodnego 
ślad tensora napięć wiąże się z ciśnieniem p poprzez związek 


/25;30/ 


dn 3 Tr (PI = okgl - ACACHCH GŁ me 1) 


x Po (%2 0000) A 


Prawa strona (3.11) nie żależy od wyboru C}. 


3.2. Reguły sum i napięcie powierzchniowe 
Rozważmy układ w którym niejednorodności gęstości wystę- 
pują tylko wzdłuż osi z. Ze względu na jednoosiową symetrię 
obrotową i symetrię zwierciadlaną względem płaszczyzn równole- 
CH 
głych dodz tensor napięć E (z) ma tylko dwie niezależne składo- 


we: transwersalną Pr (Z) i normalną PN (z), i równa sie 
P= pri) (Gre 6, + AeA) + put) bø , (3.42) 


gdzie z o: 8. są wersorami odpowiednio wzdłuż osi x,y,z. 


Y 
Równanie hydrostatyki (3.9) upraszcza się do postaci: 


PU = -| do fizo) Wego) (343) 


Całkując (3.14) od punktu z, do nieskończoności otrzymujemy 
równanie wiążące ciśnienie p=py (99) z profilem gęstości 


P= Pu(zo) -fez fdo Plz) A aj (3.14) 


Zbadajmy (3.14) w układzie twardych cząstek. przy twardej 
ścianie umieszczonej w z=z,=0. 


Potencjał zewnętrzny ściany wyraża się wzorem 


too  ZŁZm(w) (3.15) 


0 Z > Zm(w) 
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Vea 0)= 


2,25 


gdzie z, (uw) jest minimalną odległością na jaką twarda 


cząstka o orientacji C9 może zbliżyć sie do twardej ściany. 


Ponieważ funkcja 


MG) = PA(zio) EX på B Vare (iw)j (3.16) 


jest ciągła w punkcie z=z„ (W), 


głoska) Noe E) 4 2 Peder) 48C-2mtd), (3.17) 


a Py(9)=0, więc ciśnienie p w równaniu (3.14) wyraża się 
wzorem /31/ 


P= KaT |do m (Zm), (348) 


Powyższa reguła sum nie była dotąd podana. Dla twardych kul 


o średnicy O” redukuje się ona do znanego wyniku /32/: 
p=kgr m($) (2.19) 


wiążącego ciśnienie w układzie z gęstością kul na ścianie. 
W układzie z powierzchnią, wielki potencjał termodynamiczny 
2 (2.18) rozkłada się na sumę dwóch części: powierzchniowej 
Als i objętościowej L2p=-pv, 


Q=Q5-pV. (3.20) 


Jeżeli zwiększymy rozmiary liniowe układu w kierunku wersora 


A 
e leżącego w płaszczyźnie xy o czynnik dæ to zmiana Q „ALL, jest 


równa pracy wykonanej nad układem. 


AN=LAJETP BJ dz= „CA faraz, (324) 


gdzie A jest polem powierzchni. Korzystając z ekstensywności 
SL otrzymujemy że 
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Q=-AjprZ)dz , (3.22) 


Łącząc (3. 14), G. 20) iB. 22) dostajemy wzór na fs, 


zie ą5 = (dz [A2)-pre)] - [dz [do z 20 
sd Vo xe (0) (2. 2 3) 
YZ ) 


który dla twardej ściany redukuje się do następującej 
reguły sum /31/ 


Qs- faz [Ps 2)-plz)]+ kal Ja Zm (00) (Zu) 0) (3.24) 


Pierwszy wyraz we wzorze (3.23) będący wkładem do energii 
powierzchniowej pochodzącym z oddziaływań międzycząsteczkowych 


jest napięciem powierzchniowym, 


SE 32 
y= fdz CPne)- pre)]. (3.25) 
Korzystając ze wzorów (3.9)- (3.10) można zapisać p(z) i 

Pr (z) w następującej postaci /30/: 


Pu C=) = ker gen) - $ [a faros feo za EEK), 
«fax fa (z-dzia , Ta, 40), (3.26) 
Pr (2)= kalple) — 4 4 fanz jaw jaa, EE og. aS 
+ gn ale RD] ipimie (3. 27) 
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przy założeniu liniowego konturu Cj, tj 1=- drynt 1 , Ośd<4d 
Napięcie powierzchniowe y jako wielkość mierzalna nie zależy 


od wyboru C12 i jest równa /30/ 


(= 4 faz Jama [day [dd 9 (2 42 KN + 
+ (2 bi = 22a Ż,, | (2,0%, c (3.28) 


wzór (3.28) jeśt rozszerzeniem wzoru Kirkwooda-Buffa /33/ na 
przypadek cząstek anizotropowych. Aby stosować wzór (3.28) 
musimy znać 92 dla układu niejednorodnego. W praktyce jeżeli 
znamy przybliżony funkcjonał Qy (2.20), to rozwiązanie 
równania (2.21) na $Cra) pozwala nam obliczyć Qs a tym samym 
Yy Bezpośrednio ze wzoru 6.20). Ze wzorów (2.11) i (2,18) wy- 


nika po krótkim rachunku równość 


E wo 
rer ać Sila” (3.23) 


Całkując (3.29) po całej przestrzeni dostajemy związek 

2:2) eq 

BR- Afa faod, 1330 
który dodajemy stronami do równania 


By >EN, (331) 


W wyniku otrzynujemy /patrz (3.20y następującą regułę sum 


(4) A ha" - faz{ fdu tzw) — gå. (332) 


która dla twardej ściany wiąże adsorbcję , D 


= Jaz I |(dwgzw) - pf (3.33) 


z Ne zgodnie z równaniem adsorbcji Gibsa /28/. 


http://rcin.org.pl 


- 28 = 
4. TEORIE PRZYBLIŻONE. 


Typowa gęstość substancji ciekłokrystalicznej złożonej 
z cząstek o objętości v, jest rzędu voz; jest więc zbliżona 
do gęstości maksymalnego upakowania. Przy tak dużych gęstoś- 
ciach istotną rolę w określeniu struktury faz ciekłokrystalicz- 
nych odgrywa krótkozasięgowy i anizotropowy potencjał odpy- 


chający /34/. Potencjał ten możemy modelować twardym rdzeniem, 


4.1. Model potencjału oddziaływania i objętość wykluczona 
Zakładamy, że układ składa się z twardych wydłużonych 
: cząstek o symetrii cylindrycznej. Potencjał oddziaływania 


pary cząstek wyraża się wzorem 


ki przenikają się 
T-T U) ag OFr ó 4,4 
20% AI 2410.) jo -"- nie przenikają sie. 
Wówczas funkcja Mayera fo (2.14) równa się odpowiednio -1l i O. 
Onsager /35/ wprowadził objętość wykluczoną Vo (04,02), 


-Vo (wqicd,)= fdr fi Tia Wy, O>), (4. 2) 


którą oblicza się przez całkowanie po położeniach /środków mas/ 
cząstek przy ustalonych orientacjach (xy, Wy . Obszar rj. taki 

że f,40, przy ustalonych orientacjach obu cząstek, tworzy 

pewną bryłę, zwaną bryłą wykluczonej objętości, Rys 44.1), 4.2) 

i 6.3) pokazuje sferocylinder i bryłę wykluczonej objętości okreś- 


loną przez Onsagera dla pary sferocylindrów, 


pdłsjera cylinder 
L 


fp 
wz 
rys 24.1) sferqryliafistrg.pi 
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rys. (4.2) rzut bryły wykluczonej objętości, 
dla pary sferocylindrów, na płaszczyznę podstawy 


rys. (4.3) graniastosłup prostokątny o podstawie 
rombowej /część O na rys. (4.2)/ 
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Bryłą tą jest graniastosłup prostokątny o podstawie rombowej 
/rys.(4.3/ otoczony połówkami cylindrów /I-IV/ i wycinkami kul 
/1-4/ o promieniach D/2 /rys.4.2/. Kat rozwarcia rombu (a, 

jest kątem pomiędzy osiami symetrii cząstek. Za pomocą wektorów 
abueg i T można sparametryzować całą bryłę wykluczonej objętoś- 
ci metodą zaproponowaną przez A, Poniewierskiego /36/. Jeśli 


A A 
(04 UO), oznaczają wersory osi symetrii częstek to 


= 2 Lö (4.3) 
b= z Ló, (4.4) 
g =D (63,x6, V aA] (4.5) 
2 =p GQ, x (Oy A) [dy x (6 xd, )| (4.6) 


D 
{=D Ôx (AxA, )/ | Byr (63, x60)| (47) 


Mektor r zawarty wewnątrz bryły wykluczonej objętości można 


zapisać w następującej postaci: 
0: x=40 +tb +rę , Aśsjtyr <Ą (4.8) 


T: T= & +ób + t (C cosg + £ sing); AL 46541, (4.9) 
0<£4<1 


Ni T=b +óq-E(C toy + Gain) (4.40) 
M: T=-g tab -t(ęwóg + fan) (4.44) 
L=- +40, +t(Cuop+g sing) (4.42) 
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1: rsąrb+t(foino edy + TksidanY + (1,13) 
+£000Ø) , OśŁ<Ą , O<B<T, OSY<T-Q, 
2:T= -q+b + Ł(-eoinB cosy - FE Da sin0sinp+ 
+ 046), ( 4.44) 
3: T=-Q-h +t (f sinb osy - 2P bninBainy + 
ETT (4. 15) 
U: T=Qa-b + t (e ambeng + EPa qsimBoiny + 
+Cd0), (446) 


iza Onsagerem dostajeny że 
Vo (01,0, ) =8%, + ZLD|simey, | (4.17) 


gdzie v ATE + 10? jest objętością sferocylindra. Tak okreś- 
lona funkcja AC z (WW, ) zostanie wykorzystana w rozdziale 
5 å 6, 

Z rysunku (4.2) wynika, że gdy © „+50 /cząstki są równoległe/ 

to bryła wykluczonej objętości redukuje się do sferocylindra 

o długości 2L i średnicy 2D. 


Istnieje również inny sposób opisu funkcji fy. Zauważmy że 


fa (%-12,04,02)= å O(0(fa hc) ż [2 |) (418) 


A 
gdzie ©) jest funkcją skoku Heaviside'a, a O (hizi Gy 1 Wz ) 


jest odległością najmniejszego zbliżenia dwóch cząstek. 
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Aby obliczyć Olaa wy) wystarczy sparametryzować powierzchnię 

bryły wykluczonej objętości, ponieważ z definicji Olaz yuh) 
p 

jest równa odległości od środka bryły do jej powierzchni, 


4.2. Przybliżenie niskogęstościowe Onsagera i jego uogólnienie 
W 1949 roku L. Onsager wykazał, że w układzie twardych 

sferocylindrów następuje przejście fazowe I-N /35/. Praca Onsagera 

oparta była na następujących założeniach: 

/i/ cząstki oddziałują tylko potencjałem twardego rdzenia 

/ii/ układ ma małą gęstość, EYKA 

/iii/ cząstki są długie LD. 

Wyrażenie ef (0) Ac) określa gęstość cząstek o określonej 

orientacji zawartej w kącie bryłowym AW wokół w . Onsager 

potraktował więc układ jako płyn wieloskładnikowy w którym 

dany składnik odpowiada określonej orientacji cząstek. Korzysta- 

jąc z założenia /ii/ oraz z rozwinięć wirialnych dla układów 

wieloskładnikowych otrzymujemy w granicy Aw>() następującą 

postać energii swobodnej /35/: 


ka -Jiofofuutgo- + $faadaktawy 4 
fD) +08 


Pierwszy człon w równaniu (4.19) opisuje entropię rotacyjną 
wziętą z przeciwnym znakiem; jest ona największa dla rozkładu 
izotopowego 4(w)=4/yT . Z kolei drugi człon jest entropią transla- 
cyjną wziętą z przeciwnym znakiem; jest ona największa gdy 
cząstki są równoległe /patrz (4.1773 ponieważ wówczas przestrzeń 
dostępna dla cząstek jest największa. Minimalizując F wzgledem 
fi) otrzymujemy nieliniowe równanie na flw) „ które ma dwa roz- 


wiązaqia: izotropowe i anizotropowe. To równanie można rózwię- 
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zać rozwijając flo) w szereg wielomianów Legendre'a /37/: 
H= dr (4 rå av (cos0)) (4.20) 


Dla rozwiązania de otrzymujemy ał0. Przyrównując 
4 
ciśnienia, E 7 IW oraz potencjały chemiczne ,q=f + F 
N 
obu faz KG 3 gęstości fazy izotropowej i nematycznej, 
r A Pu , Oraz nematyczny parametr porządku Q we współist- 


nieniu: 


2 gu ŽD 2 (Pu-$r) 
Z --406, W26, Tap” 0264, (4.24) 
Q=03814, 


W granicy L->2æo gęstości fr. Pu dąży, do zera, wydaje się 
więc, że obcięcie rozwinięcia wirialnego na drugim wyrazie 
jest poprawne w tej granicy. Granicę Q% 30, 5900, QV 5 = const 
będziemy nazywać granicą Onsagera. Model Onsagera uogólnimy 

na przypadek układów niejednorodnych /18,19/. Z równości (2.14) 
oraz (2.36) otrzymujemy że Cozfo w przybliżeniu niskogęstoś- 
ciowym. Wstawiając to wyrażenie do wzoru (2.43) i kładąc 920 
otrzymujemy niskogęstościową postać funkcjonału energii swobod- 


nej dla układu niejednorodnego i anizotropowego: 


FTS] = [ardo pero) Un Nytt) — 4 ] 


(4.22) 
- 4 fay dw År dw, O(GW)P(5 2) fa 510003), 


która dla Pro) oflo) redukuje się do wzoru Onsagera. 
Przybliżenie Onsagera jest niewystarczające do opisu gęstych 
układów ciekłokrystalicznych. Zostało pokazane w szczególności, 


że w ramach tego przybliżenia nie otrzymujemy przejścia fazo- 
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wego N-SmA /19,24/. Niemniej model Onsagera dostarcza wskazówek 
potrzebnych do sformułowania teorii gęstych układów, Daje ścisły 


opis przejścia fazowego N-I w granicy pp =const, Ľ/D > oo, 


4.3. Teoria: SDA niejednorodnego płynu twardych kul 


Najprostsza teoria mikroskopowa, którą stosuje się w opi- 
sie niejednorodnych płynów prostych ,to model GvDW /generalized 
van der Waals/ /38/. W tej teorii Jẹ” (2.43) jest sumą wkładów 
pochodzących od długozasięgowych sił przyciągających i krótko- 
zasięgowych sił odpychających. Gęstość funkcjonału energii 
swobodnej (2.44) związanej tylko z siłami odpychającymi przy- 


bliża się równaniem 


NY (£ ;E96QI) = Apr CD) (4.23) 


gdzie AY pon jest gęstością energii swobodnej jednorodnego 
płynu twardych kul i funkcją lokalnej gęstości Q(T) e Przybli- 
żenie to jest słuszne jeśli skala przestrzennych zmian gęstości 
jest dużo większa niż skala krótkozasięgowych korelacji w ukła- 
dzie. Tak jest np. w zjawiskach zwilżania /39-48/, kondensacji 
/28,38,49,50/ i nukleacji /49,50/. Istnieje jednak wiele sytu- 
acji, kiedy powyższe założenie nie jest prawdziwe np. w układzie 
twardych kul przy twardej ścianie lub w przejściu fazowym 
ciecz-ciało stałe, a więc wtedy gdy krótkozasięgowe korelacje 
prowadzą do pojawienia się silnych niejednorodności na obszarach 
rzędu rozmiarów twardej kuli. Tego typu zjawiska są badane w 
ramach teorii SDA /51-62/. Poniżej omówimy -jedną z wersji SDA 
zaproponowaną przez Tarazonę i Evansa-/51-53/ oraz Curtina 

i Ashcrofta /55-59/. 


W SDA zakłada się, że AV /röwn.(2.44Y może być przybliżona 
przez Ahom dla układu jednorodnego o. odpowiednio jednak 


http://rcin.org.pl 


- 35 - 


dobranej, chwilowo nieznanej gęstości 90) tjs 


AV(x 80) =AY(FO). (4.24) 


AP hon jest już zwykłą funkcję Ọ (x). Efektywna gęstość pr) 
jest z kolei funkcjonałem (x) » Ponieważ dla twardych kul 
AV hon jest jednoznaczną funkcją gęstości gwarantuje to ścis- 
Łość (4.24) przy odpowiednim wyborze funkcjonalnej zależności 
ox) od o(x) tzn. UCZYC) . W praktyce musimy 
podać przybliżony sposób obliczania P(X) . Zakładamy, że 


F= [woeri pæ) endr. (4.25) 


Waga w jest funkcja g(x) i ma zasięg bezpośredniej funkcji 
korelacji c /patrz równ. (2.437, a więc szybko dąży do zera 
dla Ieho /6 -średnica twardej kuli/. Zasięg w odzwierciedla 
skalę nielokalności A (G[g(5])/patrz 8 2.4/. Funkcja w musi 


spełniać warunek normalizacji, 
far’ wie tipo), 


ponieważ w układzie jednorodnym 9=9 . Efektywna gęstość 
zmienia się wolno w przestrzeni /58/ i nie przekracza nigdy 
gęstości maksymalnego upakowania w przeciwieństwie do lokalnej 
gęstości (7) /57/. 

Aby w pełni określić >. daga musimy znać wagę w oraz ŚP pon” Obie 
funkcje możemy powiązać z c dla układu jednorodnego. AP om 
łatwo odtworzyć korzystając z relacji (2.37) dla twardych kul, 
a równanie na w otrzymamy ze związku (2.32) /55/. Łatwo spraw- 


dzić, że dla układu jednorodnego o gęstości Q otrzymujemy 
Ca (15-4 9) = -2p AKam(9) WOEN 59) -B9 Aal) fari wieć 5;9) wiri) 
I myl. 3 
- Bg ng) Jo vily-xå 9) MG HDs pig) 5 p) (4.26) 
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AV hon i AY on oznaczają pochodne względem 9 « Równanie to 
określa w(c=r";9) jako funkcję r-r' i jako funkcję  . W ukła- 
dzie niejednorodnym podstawiamy w miejsce 9 pl) i to w pełni 
definiuje w(c-r"; (7) + Własności jednorodnego płynu twardych 
kul niezbędne dla SDA są już dobrze poznane /58,61,53/ i obok 
takich wielkości jak c. można znaleźć nawet wielocząstkowe 
funkcje korelacji /64,65/. Przy badaniu układu ciekłokrystalicz- 
nego trudności pojawiają się już przy badaniu własności płynu 


jednorodnego. 


4.4. Teoria SDA niejednorodnych układów złożonych z twardych 
anizotropowych częste 


Sytuacja jest znacznie bardziej skomplikowana przy kons- 
truowaniu teorii niejednorodnych układów złożonych z twardych 
anizotropowych cząstek, ponieważ nawet własności fazy izotro- 
powej są słabo poznane. W szczególności nie znamy ani A Phon? 
ani APTN w„)dla jednorodnej fazy ciektokrysta—licznej. 
Dlatego też chociaż SDA dla twardych kul dostarcza nam wska- 
zówek, jak budować JF to jednak nie będziemy mogli skorzystać 
z przepisu na obliczanie A+ ków / z równ. (2.37 i w / z równ. 
(4.26). 


Wychodząc z równań 


Figo Falgu] + Foto], (4.23) 
Fig [gt50)]= kgl fdrdw pt) fun Ko (rio) -4$, (4.26) 


przyjmujemy wyrażenie typu SDA na F* 


FO [giro] jar gir) Meg (POD) (4.29) 
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"gdzie NY ref Jest gęstością energii swobodnej jednorodnego 
płynu anizotropowego z pewną orientacyjną funkcją rozkładu 


fref (w) „W pracy rozważymy dwa przypadki: 


1/ f,,ę=1/4, 4.30 
11/ fop Ar] ERD- +STAD+DJ. A 


Z równania (4.29) wynika, że cała zależność F” od orientacyj- 
nej funkcji rozkładu f(r,w) zawarta jest w gęstości DU). 

Z analizy modelu Onsagera /patrz komentarz do wzoru (4.19) 
wynika, że im bardziej anizotropowa jest orientacyjna funkcja 
rozkładu f(r,w) tym "słabiej" cząstka w punkcie r odczuwa 
obecność innych cząstek. Wskutek tego F*, jako część energii 
swobodnej wynikającej z oddziaływań, ulega zmniejszeniu przy 
zwiększaniu anizotropii. Fakt ten można uwzględnić poprzez 
PL) tzn. zależność P(X) ód f(r,w) musi być taka aby O CL) 
malała przy wzroście anizotropii t(r.w) . Załóżmy przez analo- 


gię do uogólnionego modelu Onsagera (4.21), że 


(5)= Jaz donde, m(v7-17,0040,) fga) pt) Hz), (4.34) 


gdzie w jest funkcją wagową. Powinna.ona spełniać warunek unor- 


Jag dead, WOT 04100.) fręlua) freęlw)s4, (4.32) 


ponieważ w granicy f(r, w) > fg (W) 9=© . Z analizy SDA 
dla twardych kul oraz ze wzoru (2.43) wiemy, że waga u pówinna 
mieć zasięg potencjału. Najprostsza funkcja spełniająca powyż- 


szy warunek oraz równość (4.32) ma następującą postać: 


) 


W(5-%10,0,)= ~ BG 5 100) (4.33) 


Ne: 
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ref 


gdzie Ba jest drugim współczynnikiem wirialnym płynu jedno- 


rodnego, 
= -$ dnddyda, folge) Srel) reg le), (4.34) 


Na koniec zaktadamy AV ret w postaci 


At (9)= Bro + (AVM)- Hm), (4.25) 
gdzie 1=$% No -objętość cząstki/; A= S kaj jest gęstością 


energii swobodnej Carnahama-Starlinga /66/. LU vor zostało 
tak dobrane aby w granicy L/D Q odtwarzać energię swobodną 
twardych kul a w granicy 9v,> O model Onsagera. Skądinąd 
wiadomo /67/ że równanie stanu Carnahama-Starlinga dość dobrze 
opisuje izotropowy układ krótkich sferocylindrów. Jak później 
pokażemy w granicy Onsagera, SDA redukuje się do modelu Onsa- 
gera, zaś w granicy L/D O odtwarzamy najprostszą wersję SDA, 
dla płynu twardych kul, w której waga w, dana wzorem /51,52/- 
lim. (5-600) = 3 © (V- I] ) (4.36) 
Lp>0 4T 
nie zależy od gęstości. Rola jaką układ jednorodny odgrywa 
w teorii SDA dla cząstek anizotropowych jest taka sama jak dla 
twardych kul. Energia swobodna związana z niejednorodnościami 
jest zawsze obliczana względem energii swobodnej układu jedno- 
rodnego. We wzorze (4.30) przypadek /ii/ pojawie się wtedy, gdy 
będziemy badali niejednorodny układ idealnie uporządkowanych 
cząstek. Rolę fazy izotropowej przejmie wtedy idealnie upo- 
rządkowana jednorodna faza nematyczna. 


Korzystając ze wzoru (2.32) policzmy cz W tej teorii: 
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- Kel Ca (14,%0% 2.) = AP re (515) - >, - Å ) + MAGA x 


S90) j "AA SUT) SPO) 
R + fargen A, (F) N gun) * 


(gy) & 2 437 
+ jar 9) Ak GUD) Fer ? (437) 


gdzie Av fe). Av! rad (26) oznaczaja pochodne czastkowe 


AY pef względem 900) .„ Wprowadźmy funkcję pomocniczą 


P (5,0%) = [dr dw, WT 04, 0) PCR 0 ) (4.38) 
i ekorzystajmy ze związku (4.31). Otrzymujemy 
SP 5 å 
ork)” (%- 0 -5(7) + $ (210)| + Śr do giro)» 

x WO), (4.39) 
So) __8G-0) Søm) SEn) SPA) + 
aaaea) SD ra) PL) 39070) 


EG 2, + 8660 vi (0-5 jua). (4.40) 


Stad 


-klea lagena) Abia (QD)MGG 10000.) + 
AV (GE) AF, ea) -2)/9(5) + 
+Z AG) APG) dw Mag Os DET 
+ far dodw POD AG PG) WCG 04,6) Waga) 
x (ra) f(r) , (4.44) 
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gdzie Aglxu)= 9050)- Pf „w ten sposób w i AY at określają 
ca. Łatwo sprawdzić, że równanie (2.37) wiążące A Weg i c 
spełnione jest tozsamo$ciowo. 

Teoria SDA przedstawiona w tym paragrafie może być stosowa- 
na do badanta układu cząstek o dowolnym AE PE Cała infor- 
macja o tym jaki układ badany zawarta jest w dwóch wielkościach: 


VE i foe 


W następnym rozdziale zbadamy zjawiska powierzchniowe 


w fazie nematycznej w granicznym przypadku modelu Onsagera. 
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5. ZJAWISKA POWIERZCHNIOWE W CIEKŁYCH KRYSZTAŁACH 
ZASTOSOWANIE UOGÓLNIONEGO MODELU ONSAGERA 


W rozdziale tym zostaną przedstawione wyniki prac /68-70/. 


5.1. Obszar. międzyfazowy faza izotropowa - faza nematyczna 
W uogólnionym modelu Onsagera (4.22) „(2.20) równanie na 


profil równowagowy (2.21) ma postać 
Fa d (51) 
PO) Ea Pix (m) tfd do falnar) (S 


Bezpośrednie rozwiązywanie równania (5.1) jest zbyt skomplikowa- 
ne i musimy odwołać się do przybliżeń, Jednocząsteczkową zredu- 
kowaną funkcję rozkładu odpowiadającą współistnieniu fazy izo- 
tropowej nematycznej przybliżamy w następujący sposób: 

olw) = ( pyr  z<O, 

PNF) 2>0, 

kładąc Vaxt=0- Jest to tzw. przybliżenie ostrej powierzchni 
/sharp interface approximation/. Lewą półprzestrzeń zajmuje 
faza izotropowa o stałej gęstości fr , a prawą faza nematyczna 
o stałej gęstości PH i orientacyjnej funkcji rozkładu HRD) 
(4.20) , przy czym zakładamy, że kierunek fi jest stały w pół- 
przestrzeni z) 0. "Direktor" ń jest teraz traktowany jako 
zmienna od której zależy wielki potencjał termodynamiczny 
NQ [elr] (4.22), (2.20) . Ponieważ żaden kierunek w płaszczyź- 
nie xy nie jest wyróżniony (Q zależy jedynie od kąta jaki â 
tworzy z normalną do powierzchni tzn. od kąta pochylenia 0 
/tilt angle/. Równowagowy kąt pochylenia of odpowiada minimum 
SE « Objętościowa część wielkiego potencjału termodynamicznego 
Q, nie zależy od n. Wynik ten jest ogólny, W naszym przypadku 
/wzór (4.197 wynika to z faktu, że objętość wykluczona AOC? 


ma pełną symetrię rotacyjną. W nieobecności pól zewnętrznych 


powierzchniowa część wielkiego pôtencjažú termodynamicznego na 
http://rcin.org.p 


a 4Ż = 


jednostkę powierzchni jest równa napięcłu powierzchniowenu, 

Y y- (3.25). Stąd otrzymujemy, że równowagowy kąt og 
odpowiada minimum napięcia powierzchniowego ý= (0-12, )/A 
/patrz(3.20/. Odejmując od wielkiego potencjału termodynamiczne- 
go QL 6.22), (2.20) część objętościową Kor (4.19), (220) 1 ko- 
rzystając z (5.2) otrzymujemy napięcie powierzchniowe faza 


nematyczna - faza izotropowa, ywr ; 


Dy, = NI = Z ż kar | faor dod, |- (5) ffin ER 
AAC RaZ ỌN Ja,(| faz V zz|©4£0>) — V (ov ',)flajfia) 
+2 On (ŻE z) ffa, V(lzzl du) two] p (5.3) 
V (12420410, )= - (drz fa (14210400) . (5.4) 


Piz jest rzutem r.» na płaszczyznę prostopadłą do osi z. 


Meżygykie całki występujące w (5.3) przekształcamy do postaci: 
ji [dza Vbolowa)=faz, zaVlkelcąch)=V, (c); (55) 
sted otrzymujemy 

Sn C) = pi (A) + (A) + Nå > (5:6) 
gdzie 


ga(R)/keT= GOL (deg de, Vi (09409, JE lea) , (52) 
ve) IkęT - SÅ dody dcą V (du) FP), (58 
a/kar = — 4 (LEN [dey doa V (udo), (53) 
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Zarówno ya jak i y2 zależą od ñ, ponieważ V,(04, Wa) /dodatek AJ 
nie ma pełnej symetrii rotacyjnej. Fizyczna interpretacja 
MADE jest następująca: 

Cząsteczka ciekłego kryształu należęca do fazy nematycznej 
/lub izotropowej/ i znajdująca się w pobliżu powierzchni roz- 
działu faz z=0 ma mniej najbliższych sąsiadów należących do 

tej samej fazy w porównaniu z cząsteczką w objętości. Wzrost 
entropii translacyjnej związany z tym efektem opisują odpowied 
nio wyrazy - y2 i - p- Ujemne wkłady Va i 65 do napięcia 
powierzchniowego fer skompensowane są przez dodatni wkład p 
pochodzący z' bezpośrednich oddziaływań cząsteczek należących 
do fazy izotropowej z cząsteczkami należącymi do fazy nema- 
tycznej. Napięcie powierzchniowe Yi-r musi być dodatnie 

aby układ z powierzchnią był stabilny. 

Obliczenia numeryczne zostały wykonane dla układu sferocylin- 
drów o długości L/D=5,8,10,15,20,25,100. Gęstości PN i fr 
współistniejących faz (4.21) oraz fl) (4.20) zostały wzięte 

z pracy Lashera /37/. Całki wielokrotne policzone zostały me- 
todą Monte Carlo. 

Na wykresie (5.1) przedstawiamy y2 w funkcji Q. Funkcja %2 

ma minimum dla 0,= 0” niezależnie od długości L/D /röwniez 

dla L/D=15,25,100/. Z kolei Vi /wykres 6.2)/ma minimum dla. 

64, =900. Wyniki te można wytłumaczyć w następujący sposób: 
Jeżeli tylko połowa przestrzeni jest wypełniona fazą nematyczną 
a druga połowa jest pusta to wydłużona czasteczka, której 
środek masy znajduje się na powierzchni z=0 ma tym większą 
swobodę ruchu w porównaniu z cząsteczką w objętości im bardziej 


wystaje na drugą strong. Dlatego %2 ma minimum dla ustawienia 
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prostopadłego cząsteczek względem powierzchni. Jeżeli pustą 
półprzestrzeń wypełnimy teraz fazą izotropową to cząsteczka 

fazy nematycznej na powierzchni zyska pewną liczbę sęsiadów 

z fazy izotropowej w porównaniu z cząsteczką w objętości. 
Zmniejszenie entropii translacyjnej tej cząsteczki /w. porówna- 
niu z cząsteczką w objętości/, związane z tym efektem i opisy- 
wane przez -Y4 , jest tym większe im bardziej cząsteczka wystaje 
na drugą stronę, Stąd ya ma minimum dla ustawienia równoległego 


cząsteczek względem powierzchni. 


.09 309 % 60° 909 


-04 


nD? 
kgl 


rys.6.1) aw funkcji kąta pochylenia 
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rys.6.2) Vi w funkcji kąta pochylenia 
Przeciwstawne tendencj /rys .B.1) i (5.2) prowadzą do uporządkowa- 
nia cząstek na powierzchni pod kątem Op 7=607+40 dla którego 
YnI= Mejer å ma minimum /yg 6.30 dla wszystkich L/D. Napiecie 
powierzchniowe jest dodatnie i układ jest stabilny. 
Na wykresie 6.4) przedstawione jest równowagowe napięcie powierz- 
chniowe ACE funkcji L/D. Wielkość ta malejedo zera w granicy 
Onsagera /54.2/ ponieważ VoQN „Wór: dążą do zera w tej granicy. 

Gęstość fazy nematycznej ON jest wielokrotnie wieksza niz 

gęstość pary fu, wiec możemy uznać (2 za główny wkład dò napię- 
cia powierzchniowego faza nematyczna - para pochodzący z oddzia- 
ływań twardych /wzór (5.6)- (5.9) z Qr zastępionym przez /. 
Oczywiście w układzie cząstek twardych nie następuje zjawisko 


kondensacji, niemniej poniższe wyniki dają jakościowe wskazówki 


co do wkładu twardych oddziaływań w zjawiska zachodzące na 
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rys .6.3) napięcie powierzchniowe w funkcji kata 
pochylenia 


002 


004 0.1 a2 


Au 


rys (5.4) równowagowe napięcie powierzchniowe 
w fuhkcji D/L 
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powierzchni ciekły kryształ - para /powierzchnia swobodna/. 
Na powierzchni swobodnej twarde oddziaływania uprzywilejowują 
ustawienie prostopadłe / 9, =0°, rys 6.1). W punkcie potrójnym 
/współistnienie fazy nematycznej, fazy izotropowej i pary/ 
napięcie powierzchniowe zmienia się skakowo od napięcia po- 
wierzchniowego faza izotropowa -= para, fry= 53 . do napięcia 
powierzchniowego iza nematyczna = para, Java: 

Otrzymujeny 


Ayw = Ya (0) -ys <0 (5.40) 


Na zakończenie warto zauważyć, że wzory (5.9-6.9) można również 
otrzymać stosując bezpośrednio wzór Kirkwooda=Buffa (3.28), 
przy założeniu niskogęstościowej postaci dwucząsteczkowej zre- 


dukowanej funkcji rozkładu 
POK.) = 905 nd) PE) (EM) 


gdzie (y W) dane jest wzorem (5.2). 


5.2. Obszar międzyfazowy nematyk - twarda ściana 


5.2.1. Powierzchniowa część wielkiego potencjału termodynamicz- 
nego 


Jako pierwsze przybliżenie równowagowego profilu (70) 
/równ.(5.Jy dla nematyka ograniczonego płaską ścianą umieszczoną 


w z=0 przyjmujemy HW 
Yw) = ppf (w) EP Ver (700) (542) 


gdzie Peflu) jest jednocząsteczkowę zredukowang funkcją roz- 
kładu w objętości, a potencjał twardej ściany VE. dany jest 
wzorem (3.15). Przybliżenie (5.12) można otrzymać przy itera- 


cyjnym rozwiązywaniu całkowego równania (5.1) ; rozpoczynając 
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iteracje od rozwiązania objętościowego, Qof). otrzymujemy 
w pierwszym kroku (5.12). Dla sferocylindrów odległość naj- 


mniejszego zbliżenia występujące w (3.15) równa jest 
4 
Zm(w)= 7 (L|c00|+ D), (5.43) 


gdzie 0 jest kątem pomiędzy osią symetrii sferocylindra a nor- 
malng do powierzchni, 

Powierzchniową część wielkiego potencjału termodynamicz- * 
nego Qs obliczamy korzystając ze wzoru (3.20) tak jak w § 5.1. 


Po prostych przekształceniach otrzymujemy D na jednostkę po- 


wierzchni A: 


Jk" EE, (5:44) 


gdzie S jest entropią powierzchniową, © jest adsorbcje 
/wzör 3.33/ równą 


l= - og deo Zmed FC). (5,45) 


Entropia powierzchniowa składa się z dwóch członków: 


entropii rotacyjnej i translacyjnej, 


3 Ów * 63 (5:46) 
gdzie 
Śrę/kg = |do zit) Ln[4nf(0)] (co), (5,47) 
Jer = = $> + tr ) (518) 


e =-9 Jao f0) zmt) fA- dn Cpa] f, (549) 
SE. $% $ [dada Flor) Fl) zmia) Voley) + (5.20) 


+% TATAN - Zu) zywy) Vane) 
Zmol)-Zml h) 
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Funkcja v( lz42| Wy) Zwzór (5.47 jest obliczona w dodatku B. 
Entropię rotacyjna, Srot: zdefiniowaliśmy tak aby znikała dla 
fazy izotropowej. sid jest zawsze ujemne ponieważ twarda ściana 
ogranicza przestrzeń dostępną dla cząstek. Również [* jest ujem- 
ne, ponieważ liczba cząstek przy ścianie ulega zmniejszeniu, 

ŚL jest poprawką do S,, wynikającę z oddziaływań pomiędzy 


cząstkami. Twarda ściana ogranicza liczbę oddziaływujących par 


x 
r 


również swobodę rotacji. Porównując wzory (3.24)-8.25) otrzymu- 


cząstek i dlatego st jest dodatnia. Twarda ściana ogranicza 


jemy napięcie powierzchniowe 
V= ęż - %jdoznia) fl). (5:24) 


W przeciwieństwie do napięcia powierzchniowego faza izotropowa 

— faza nematyczna, napięcie powierzchniowe faza nematyczna = 
twarda ściana nie musi być dodatnie aby układ był stabilny. 
Dodatniość całego {ls gwarantuje stabilność układu. Obliczenia 
0508" wfunkcji kąta pochylenia dla fazy nematycznej w kontakcie 
z twardą ścianą zostały przeprowadzone dla warunków współist- 
nienia N-I (4.29-4.21) e Wszystkie całki wielokrotne występujące 
we wzorach 6.19-5.20) zostały obliczone metodą Monte Carlo 


przy użyciu programu Vegas /71-72/. 


-Spgt å -Sęp mają minima dla ustawienia prostopadłego. Z kolei 
„sid oraz -aP są minimalne dla ustawienia wównoległego. 


Efekt ograniczonej swobody ruchu translacyjnego przy ścianie 
jest dominujący /rys.5.5/. Na rysunku (5.5) przedstawiona jest 
zależność ov od kąta O, „. Dla wszystkich L/D=5,8,10,20,100 
a ma minimum w 0; =90" /ustawienie równoległe cząstek do 
ściany/. 


http://rcin.org.pl 


-= 50 = 


2 


09 309 60° 909 
t 
rys 5.5) powierzchniowa część wielkiego potencjału 


termodynamicznego układu faza nematyczna - twarda 
ściana w funkcji kąta pochylenia 


Wstawiając £0)=4[4T w równaniach (5.1e)-5.20) otrzymujemy 
powierzchniową część wielkiego potencjału termodynamicznego 
układu faza izotropowa-twarda ściana, leż; + W warunkach 
współistnienia następująca nierówność jest spełniona w naszej 


teorii: 


QH (90) /A + Yar (Og) 05 /A (522) 


dla dowolnego O, / yu- wykres 5.3/. W warunkach współist- 
nienia N-I ściana uprzywilejowuje fazę nematyczną dla Q, =90°. 
Nierówność (5.22) oznacza także, że w układzie twardych stana 
cylindrów przy twardej ścianie możemy się spodziewać całkowitegi 
zwilżania ściany przez fazę nematyczną uporządkowaną równolegle 


do ściany. Co prawda we współistnieniu N-I w warunkach zwilża- 


nia ściany przez fazę nematyczną powinna być ściśle spełniona 


równość http://rcin.org.pl 
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Q$/A + Yur =(15/A 


a to że zamiast niej otrzymujemy nierówność wynika z poczynio- 
nych przybliżeń. W przypadku ustawienia prostopadłego /które 


nie odpowiada minimum Oy otrzymujemy przeciwną nierówność, 
NW IW 
OS (OA > Yw-r (%) + QZ"/A , (5.23) 
spełnioną dla dowolnego OŁ z 


5.2.2. Gęstość i profile parametrów porządku 


W przybliżeniu (5.12) grubość obszaru międzyfazowego 1 jest 


hw 
ext 


jak i parametry porządku Q(z) i P(z) , (1.3-(.4) przyjmują wartość 
objętościową. Dla z <D/2 f(z) =0 i Q(z) oraz P(z) nie są 
określone.» Profil gęstości 2)=glzw)dw wyraża się wzorem 


On 
Qz)=< gab), rpdorjeosna gt), (524) 


gdzie O „(z) = arcos wa Az=z-D/2. 


Tensor parametru uporządkowania jest wą przez wyrażenie 


RORE 4 GL 24 (Broz; Sehr) 
Uj =< Goraj- Nig Pa Sk - a (5.25) 


gdzie (0)= 200) jest orientacyjną funkcją rozkładu /§ 2.1/. 
z 


równa zasięgowi V tj 1=(L+D)/2. Dla z>U zarówno gęstość 


EEE SEE Haiocier Aaaa 
9(z) = Qp 242, (5:26) 
Qz (2) „4 (zde 4 E (527) 

Qxx (2) = Qy(2)=-4Q,, (2), (5:28) 
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dla 04 AzśL/2 Q,,(z) jest ujemne, co oznacza, że cząstki 
układają się równolegle do ściany. Niemniej symetria w płasz- 
czyźnie xy nie jest złamana i obszar przyścienny nie ma 
struktury fazy nematycznej. Gdy faza nematyczna jest w kontak- 
cie z twardą ścianę wszystkie profile zależą od orientacji n. 
Dla uporządkowania równoległego / Q, =30%/ wzdłuż osi x otrzy- 


nujemy 


Q(z)= Oxx (2)= $ < Żain28 SP-A mu (529) 
P(z)= 5 [Qo e) -Qy 2)] -8 costa LA y (530) 


Dla Åz=0 dostajemy związek pomiędzy głównym parametrem po- 
rządku na ścianie, Q,=0(2=0/2), i dwuosiowością, Po= P(z=0/2), 


w przybliżeniu (5.12) 

B--4 (1-0) 1534) 
Gęstość plz) oraz parametry porządku Q(z) i P(z) sę przedstawione 
na rys, (5.6) - 6.2) „ Gęstość przy ścianie jest zmniejszona i 
dąży do zera gdy z => D/2. Główny parametr porządku jest większy 
na ścianie niż w objętości. Ze względu na dużą wartość Q(z) 


dwuosiowość P(z) jest mała w całym obszarze przyściennym, 


40 
o[z) 
% 
05 
0,0 O, 1,0 
242 s 
L 


rys.6.6) gęstość fazy nematycznej przy ścianie 
w warunkach współistnienia N-I 


http://rcin.org.pl 


- 53 -= 


99 
G (2) 
0.8 


0.7 


0.0 0.5 1.0 


24z 
L 


rys (5.7) główny paranetr porządku nenatyc nago 
przy ścianie w warunkach współistnienia N-I 


2az 
00 a5 L 10 

P(z) 

- 0.025 

-0.050 


rys.(5.8) dwuosiowość przy ścianie vi warunkach 
współistnienia N-I 


Konsekwencją przybliżenia (5.12) jest fakt, że profile P(z) „2(z) 
i P(z) jako funkcje 2Az/L nie zależą od L/D. 
nak hv 2(7) dla uke Ve 


Qzz (2) = faea R (0) f(0)/ 14 Jooh fu) (5. 32) 


W pobliżu | Om OWE 242 a stąd Q,,(2) ŚO na ścianie. 
Qz (z) maleje od dodatniej wartości objętościowej do Q,,=-1/2 
u Az=0. Wynika stąd, że cząstki ciekłego kryształu porządkują 
się tuż przy ścianie tak jak w fazie izotropowej (5.24-5.28), 


co jest w zgodzie z nierównościę (5.23). 
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W przybliżeniu (5.12) uwzględniamy tylko bezpośredni wpływ 
oddziaływania cząstka-ściana na profile gęstości i parametrów 
porządku, W szczególności reguła sum (3.18) nie jest spełniona, 
wstawiając (5.12) do (3.18) otrzymujemy 

P= Pb kgl. (5:33) 
Twarda ściana oddziałująęc z cząstkani znajdującymi się przy 
powierzchni pośrednio działa na cząstki umieszczone w objętości, 
ma więc nielokalny wpływ na profile gęstości i parametrów 
uporządkowania. ¥/ naszym modelu możemy uwzględnić efekty nielo- 
kalne w sposób przybliżony oraz ocenić ich wpływ na gęstość i 
profile parametrów porządku /31/. W tym celu wyraźny et» 
wystepujacy w röwnaniu (5.1) przez Ps flo). Kładąc Vaxe=0 W 


(5.1) otrzymujemy 


pu , 
a = of (4) exp (-jdradoo fa o )94 FT. (5.34) 
Wstawny (5.34) do równania (5.1) i zlinearyzujmy prawą stronę 


równania względem (260) - flo). 


Otrzymujemy równanie 
hw 
pizon) = $eflwi) eP Vere a oa) A + (5.35) 
+ far, do, fa (zaw) (Plzaeo)-Prflui)) | 


Rozwiążmy to równanie iteracyjne startując z rozwiązania 
$(z,0)= Peflu). Po pierwszej iteracji dostajemy (5.12). Druga 


iteracja prowadzi do następującego wyniku: 


Olza) = Goto) BE AW JĄ (5.36) 
+ faz, do V (faa 00) elor) (PD pf, 
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Drugi wyraz vie wzorze (5.36) jest pierwszą nielokalną poprawką 
do (5.12). Wyraz ten sprzęga ścianę z cząstką 2 poprzez cząstkę 
1. Przybliżenie(5.36) ściśle spełnia regułę sum (3.18) ; wstawia- 
dac (5.36) do (3.18) otrzymujemy ciśnienie w niskogęstościowym 


modelu Onsagera 
P=PskeT + 7 få [dayde Value) $$). (537) 


Na rysunkach (5.9) - 6. Eo przedstawione są profile (z). 

Q(z) i P(z) dla L/D=5, GQ; =90° . Dla porównania linią przerywaną 
zostały zaznaczone profile gz) .2(2), P(z) dane przez przybli- 
żenie lokalne (5.12). Z wykresów wynika, że nielokalny wpływ 
ściany prowadzi do oscylacji Q(z), Q (z) i P(z). Zmniejsza 

także parametr porządku Q, na ścianie i zwiększa dwuosiowość wk 
P, W porównaniu z profilem lokalnym. Płytkie minim Q(z) A OŚW 
są wyraźnie przesunięte względem maksimum p (z). Zmiana znaku 
P(z) dla 2Az/L=0.9 oznacza, że kierunek uporządkowania w płasz- 
czyźnie yz prostopadłej do kierunku "direktora" /ośx/ zmienia 
się. Wszystkie profile wyrażone w funkcji 24 z/L słabo zależą 

od L/D. 

Na zakończenie zbadajmy (5.12) i (5.36) w granicy L/D>0. /twar- 


de kule/. Otrzymujemy 


SO ED: 
sa] p, ,2>0/2, (538 


dla (5.12) oraz 


( s KER, 
PA: (pf (2-2) FiR-20] 13 8), 2 $z 
% a Z?ŹD, 


JA 


dla (5.36). http://rcin.org.pl 


16 2(z-B)/L 


Q 
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o 

S 

Q 

I o OD 

9 Aje = S S 
m 


rys. (5.9) gęstość przy ścianie, L/D=5, 04 2909 
Linia ciągia - przybliżenie nielokalne (5.36) 
Linia przerywana =- " - lokalne (5.12) 
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2(z-F)/L 


1.0 45 


0.5 


O 
3 8 8 go 
o oS jo) |) 


Q (z) 


rys 6.10 główny parametr porządku przy ścianie 
/patrz (5.9) 
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rys. (5.11) dwuosiowość przy ścianie /patrz 6.9)/ 
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5.3 Idealnie uporządkowana faza nematyczna przy twardej ścianie 


Zakładamy, że układ został umieszczony w polu zewnętrznym 
vB (w) . Korzystając z przybliżenia (5.12) łatwo pokazać, że 
powierzchniowa część wielkiego potencjału termodynamicznego wy- 


raża gle wzorem 


QR = '- 9, [do Zm (w) TA (wlw), (5 40) 


NE 
gdzie a jest dane wzorami (5.1446.20). 

b A 
Załóżmy, że Voxt=-hP (ÂA). W granicy h>oo /patrz (5.1) z Vaxt= 


b 
Vaxt/ Otrzymujemy 


Hofa) Å (6(86-4) + 6(RÓ+A1)) (544) 


Jest to przypadek idealnego uporządkowania wzdłuż kierunku n. 
W wyrażeniu (5.40) występują teraz dwa człony rozbieżne: Inf (w) 
oraz s które jednak znoszą się ze względu na (5.1). 
Korzystając z tego możemy przepisać (5.40) w postaci 


A= Zm (04) + 2 på Zml&) Vo(0£,9%) — 
79% Vi (94,4). (5.42) 


Jest to wzór na Q dla układu idealnie uporządkowanych cząstek 
pod kątem Q, do normalnej. Kat a, jest wyznaczony przez kierunek 


pola vb . Zastosyjny (5.42) do układu twardych cylindrów. Odle- 


ext 
głość najmniejszego zbliżenia do ściany jest równa, 


Zm (6) = 7 ( LIc»s6,] +DaimQ, ) (5.43) 


Bryłą wykluczonej objętości jest cylinder długości 2L i szero- 
kości 2D; V o=2 70? L. Taką bryłę łatwo jest sparametryzować i 
policzyć Vi ete (5.5) w analogiczny sposób jak to było zrobione 
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v dodatku A. Po wykonaniu wszystkich obliczeń otrzymujemy 


następujące wyrażenie /dla twardych cylindrów/ 


A (x050, + X xsin) + LES T (xog, + xine, )- 
EEE qay otga (50) 
ar Jå G (x080, + Xin, )+ TR BS (xeos, + sing) - 


Akt - Asa a | bag T X tga, -A MIANA I 
t Ż êta? Q Vega +xtgta, arcsin (RY, 


gd 1/x < 70, 


gdzie 4 =21%0 © a x=D/L. 
A została również obliczona numerycznie dla układu twardych 
aferocylindrów. Rysunki (5.12)1 6.13) przedstawiają My dla 


układu sferocylindrów i cylindrów dla Å =10.6 i 1/x=5,6,8. 


b 
ext 


/rys.5.5/ i z polem /rys 6.12) jest zachowanie Ms w pobliżu 


Istotnę jakościową różnicą pomiędzy układem bez pola V 


OŁ =90°. W pierwszym przypadku 30. lars =0 /rys „6.57 w drugim 
pochodna jest niezerowa /rys.5.l2/. Dla układu cylindrów 


20 


AD? 
AkgT 


10 


09 309 609 909 
rys.5.12 fs KERR Merocylindrów 


261 = 


AD? 
ÅkgT 
15 
10 
5=8 
05 
jo jo 90? 
09 30 8, 60 


rys. (6.12) Qs dla układu cylindrów 


rys.5.13/ o długości L/D=5 Qs ma maksimum w O. © 30% 
/ry / g 5 + 


Gdy L/D>00 maksimum przesuwa się w stronę Q, = 0°, ponieważ 
w tej granicy /Ą =const/ układ cylindrów zachowuje się tak jak 
układ sferocylindröw. 


5.4. Porównanie z doświadczeniami i innymi teoriami 
Wyniki rozdziału, dotyczące kąta pochylenia er „ Są zgodne 
z wynikami eksperymentów. Dla nCB n=5,6,7,8 zmierzona wartość 
eq ; G-50 peor 9 
O; na powierzchni N-I leży w przedziale 50 -70 /73-75/ „BŁ =60 , 
/indeks teor odnosi się do wielkości teoretycznej otrzymanej w 
g o teor o 
modelu Onsagera/, a na powierzchni swobodnej +0 , 8, =0 
/76-78/. Również dla innych substancji /MBBA, EBBA/ kat 6 na 
powierzchni swobodnej jest równy 0? dla temperatur, T, bliskich 
temperaturze przejścia fazowego N-I, Ty.» jednak w tym przypadku 
og zmienia się z T i rośnie od 0° do 150 dla TLK Ty; /79-81/. 
W naszym modelu /twarde cząstki/ zależność od T nie występuje, 


Bardzo wiele substancji ciekłokrystalicznych należęcych do grup: 
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bifenyli-, estrów-, PCH- itd. zostało zbadanych przy płaskich 
powierzchniach nieorganicznych /82/. Dla wszystkich tych 
substancji OĘ-900, ©Ę 90°. 

W doświadczeniach /76/ zbadano skokową zmianę napięcia 


powierzchniowego ciekły kryształ - para w temperaturze Ty_q 


Wynik otrzymany dla 5CB i 8CB potwierdza nierówność (5.10). 
Napięcie powierzchniowe faza nematyczna - twarda ściana zostało 
zmierzone w doświadczeniu Yokoyamy i van Spranga /83/. Jakościo- 
wa zależność Aly od 8, uzyskana przez tych autorów zgadza 

się z naszym wynikiem /patrz wykres (5.5)/. 

w licznych eksprymentach przeprowadzanych dla MBBA w kon- 
takcie z ciałem stałym otrzymywano Q =0? a nie OL =907 /84-85/. 
Jak zostało jednak pokazane w znaczącym eksperymencie Ohgawary 
i innych /82/ wynik ten był spowodowany zanieczyszczeniami próbki 
powstałymi wskutek hydrolizy MBBA. Wskutek tej reakcji powstaje 
pBA, substancja polarna, powierzchniowo czynna, która zmienia 
orientalis Vere) powierzchni z równoległej na prostopadłą. 
Również przy badaniu powierzchni swobodnej uzyskuje się sprzecz- 
ne wyniki /76/. Pomiar napięcia powierzchniowego jest długo- 
trwały i w czasie pomiaru może następować akumulacja zanieczysz- 
czeń na powierzchni. Nawet nieduża ilość zanieczyszczeń powierzch- 
niowo czynnych może dramatycznie zmienić napięcie powierzchnio- 
we. Wybór 5CB i 8CB w eksperymencie Gannona i Fabera /75/ był 
podyktowany dużą stabilnością chemiczną tych substancji. Nie 
odnosi się to do wspomnianej wyżej substancji MBBA. W wielu 
eksperymentach wymusza się określoną orientację fazy nematycznej 
przez odpowiednie przygotowanie powierzchni ciała stałego. 


Również wyniki tych eksperymentów są często niejednoznaczne 


/86-87/. Dla przykładu w metodzie pocierania /rubbing/ to jaka 
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orientacja jest indukowana na poweirzchni silnie zależy od 
substancji z jakiej składa się "ściereczka" służąca do pocie- 
rania powierzchni ciała stałego. Takich efektów nie uwzględniany 
w naszych obliczeniach. 

Wpływ silnego pola objętościowego Ve na napiecie powierz- 
chniowe nematyk twarda ściana /patrz § 5.3/ nie był dotąd zba- 
dany. Wiadomo jednak /83/, że pole elektryczne o wartości prze- 
kraczającej wartość progową dla przejścia Frederiksa /2/ jest 
za słabe aby wywołać obserwowalną zmianę pochodnej Er z war- 
tości O na wartość niezerowa. "JAP 

Powstało wiele różnych teorii mikroskopowych mających na 
celu wyjaśnienie problenu orientacji ciekłego kryształu przy 
powierzchniach. Parsons wykorzystał wzór Kirkwooda-Buffa dla 
powierzchni swobodnej /wzór (3.28) wraz z przybliżeniem Fowlera- 
Kirkwooda-Buffa /88-89/ w celu obliczenia ogå » Dominującą role 
w tym modelu odgrywa przyciągający potencjał van der Vaalsa 
zależny od orientacji. Wynik tego modelu Og 900 nie zgadza 
się z wynikiem doświadczeń, W uogólnionym modelu van der Vaalsa 
Telo da Gamy /90/, cząsteczki ciekłego kryształu są twardyni 
kulani o średnicy D, które oddziałują między sobą potencjałem 
przyciągającym. 

-A (Dr) - B (DIrÉ P (9,8) , v>D i 


5.45) 
0 ] yD 


JA (ci dy 10 ) = 


W tym modelu napięcie powierzchniowe N-I a także napięcie 
powierzchniowe N-V nie zależy od 8+ + Wynika to z braku sprzę- 


żenia zmiennych orientacyjnych i translacyjnych w (5.45), 
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Dodanie wyrazu c, ( P (©) rip ) /r"),sprzęgającego ta 
zmienne, do V,,, poprawia model /91/. Kąt 0 zależy od znaku 
C, i przyjmuje wartość 0° dla C,> O lub 90? dla c< O na po- 
wierzchni swobodnej. Powierzchnia N-I nie była badana w tym mo- 


delu. Powierzchnia swobodna została również zbada a n e 
na /92/. Uwzględniłj oni zarówno potencjał przyciągający jak i 
anizotropową część potencjału odpychającego w rachunku pertur- 
bacyjnym dla płynu twardych kul /25/. Wynik wh pracy jest 
zgodny z naszym tj. oto. Również Kimura i Nakano uwzględnia- 
ją obok potencjału przycięgającego,anizotropowy potencjał od- 
pychający /93/, który w ich modelu także uprzywilejowuje 6 209, 
Cały model zależy od kilku swobodnych parametrów i w praktyce 
przez odpowiednie dobranie tych parametrów otrzymuje się dowolny 
kąt Oz przedziału (09 ,e0%) na powierzchni swobodnej i na 
powierzchni N-I. Croxton /94/ wyprowadził wzór na napięcie 
powierzchniowe nematyk-para, ynv „korzystając ze wzoru Kirkwooda- 
Buffa i przybliżenia na Qa (mia 1092) , w którym pominął korela- 
cje kątowe. Formalne rachunki przeprowadził dla Og =0? i 

O, =90°. Z dokładnością do liniowych wyrazów w parametrze 
porządku i kwadratowych wyrazów w gradientach gęstości sprowadził 
otrzymany wzór do wyrażenie typu Cahna-Hillarda /28/. Nie prze- 
dyskutował jednak jaka orientacja jest uprzywilejowana w tym 
modelu, poprzestając tylko na podaniu formalnych wzorów na 

Ś nov (09 i y nv (20% . Swój drugi model Croxton /94/ oparł na 
teorii Maiera-Saupe /95/ i wzorze (3.20) . Również i w tym modelu 
nie przedyskutował zależności Y n-v od 8, . Zbadał natomiast 
jakościowo zależność ýy N-y 9d temperatury w pobliżu punktu po- 


trójnego otrzymując dobrą zgodność z doświadczeniem, 


http://rcin.org.pl 


-= 65 — 


Z trzech modeli: Telo da Gamy /96/, Okano /97/, Bernasconi 
i innych /98/, tylko w modelu Okano uwzglednia sig role krötko- 
zasięgowych sił odpychających w określeniu or na powierzchni 
ciała stałego. Rozważył on sta /85.2/ dla sferocylindrów w 
przybliżeniu idealnego uporzędkowania i otrzymał 6-907, 
W modelu Telo da Gamy OL zależy tylko od anizotropowego po- 
tencjażu przyciągającego ściany i może przyjmować dowolne 
wartości z przedziału (o?,907). Bernasconi i inni /98/ zbadali 
wkład do napięcia powierzchniowego wynikający tylko z anizo- 
tropowych sił van der Vaalsa i pokazali, że Og zależy od 
współczynnika załamania światła, m, ciała ścażago; Dla dużych 
m ustawienie fazy nematycznej jest prostopadłe a dla małych m = 
równoległe do powierzchni ciała stałego. Jest to rezultat 
sprzeczny z wynikami doświadczeń /82/. 

Wyniki tego rozdziału wskazują, że w pierwszym rzędzie 
uprzywilejowana orientacja fazy nematycznej na powierzchni 


zależy od anizotropowych krótkozasięgowych sił odpychających. 
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6. PRZEJSCIA FAZOWE W CIEKŁYCH KRYSZTAŁACH, 
ZASTOSOWANIE SDA. 


6.1. Krótki przegląd istniejących teorii mikroskopowych 


Istnieją dwie duże klasy ciągłych modeli mikroskopowych 
stosowanych w opisie przejść fazowych w ciekłych kryształach. 
Sa to a/ modele typu Onsagera /35/, b/ modele typu Maiera- 
Saupe /95/. W modelach z grupy a/ uwzględnia się jedynie krótko- 
zasięgowe anizotropowe siły odpychające. Z kolei w modelach 
z grupy b/ zakłada się, że główną rolę w tworzeniu struktur 
ciekżokrystalicznych odgrywa długozasięgowy efektywny /typu 
pola średniego/ potencjał przycięgający o symetrii badanej 
fazy. W szczególności pomija się anizotropię kształtu częstecz- 
ki. SDA przedstawiona w rozdziale 4 należy do modeli z grupy a/ 
i dlatego na nich skupimy naszą uwagę. Opis modeli typu Maiera- 
Saupe można znaleźć w pracach /95,99-103/. 

Pierwotna metoda rozwinięć wirialnych Onsagera była stosowana 

w pracach /104-106/ do badania przejścia fazowego I-N w ukła- 
dach twardych cząstek o różnym kształcie. i! modelu Onsagera 
zbadano także mieszaniny /107-109/, a uogólniony model Onsagera 
/8 4.2/ z potencjałem przyciągającym zastosowano w badaniu 
stabilności fazy izotropowej względem zaburzeń smektycznych 
/19,20/. Idea Onsagera została rozszerzona na przypadek 
jednorodnych układów gęstych za pomocą teorii cząstki skalowanej 
/scaled particle theory > SPT/ /110-112/. Zastosowanie SPT do 


układu wydłużonych częstek prowadzi do następującego wzoru 


na JF: 


F 9(9% dod Giwi), | (6.4) 
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gdzie g jest pewną funkcją, zależną od kształtu badanych 
cząstek, Modele te zostały rozszerzone na układy cząstek 
wydłużonych oddziałujących również potencjałem przyciągającym 
/113-116/, który uwzględnia się w przybliżeniu pola średniego. 
Również dla układów gęstych została opracowana metoda rozwi- 
nięć JJ w wielkości y= P/U- gy) /y-expansion/. Vspółczynniki 
tego rozwinięcia otrzymuje się przez porównanie z pierwotnyn 
rozwinięciem wirialnym /117/. 

Do początku lat 80-tych nie zostały wykonane żadne 
symulacje komputerowe w których zaobserwowano by przejście 
fazowe N-I w układzie twardych cząstek. Wczesne symulacje 
Vieillarda-Barona /118/ dla układu sferocylindrów nie przyniosły 
oczekiwanych wyników ze względu na bardzo powolne osiąganie 
równowagi przy dużych gęstościach. Co więcej pewne wyniki teore- 
tyczne wskazywały /19,24/, że faza smektyczna nie może tworzyć 
się w układach cząstek twardych. Spowodowało to spadek zain- 
teresowania układani cząstek twardych jako modelami ciekłego 
kryształu. 

Przełom nastąpił pod wpływem niedawnych symulacji kompu- 
terowych Frenkla, Muldera i McTague /119-120/ dla układu 
elipsoid obrotowych, Stroobanta, Lekkerkerker'a i Frenkla 
/121-123/ dla układu sferocylindrów idealnie uporządkowanych 
orientacyjnie oraz Frenkla /124-125/ dla układu sferocylindrów 
o pełnej swobodzie rotacji /nieidealne uporządkowanie/. VW pra- 
cach tych obok fazy izotropowej i ciała stałego wykryto rów- 
nież trzy podstawowe fazy ciekłokrystaliczne: nematyczną, smek- 
tyczną i kolumnową i ustalono /ilościowo/ diagramy fazowe. 

W ostatnich trzech latach powstało wiele modeli opisujących 


przejście fazowe N-I w układzie twardych elipsoid obrotowych 
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ex r 
/126-130/. Pokrótce omówimy metody obliczania F ZaProRong 
wane w tych modelach. U. Singh i Y. Singh /128/ rozwinęli F R 


iso 
2 


korelacji 0-Z fazy izotropowej założyli przybliżenie Pynna ~= 
Mulfa /131/ 


050 HS Mal 
są z ( (fu 140, w) ) (6. 2) 


funkcjonalnie wokół stanu płynu izotropowego. Dla c - funkcji 


gdzie ch S jest funkcją 0-Z twardych kul, a O - odległością 
najmniejszego zbliżenia dwóch elipsoid, (I celu obliczenia 

(og autorzy wykorzystali przybliżenie Bernc-Pechuckasa /132/. 
Baus i inni /127/ zastosowali ścisły wzór (2.43), przy 


czym założyli, że 


law) G! (I Tal; 9) dla fazy izotropowej 
Ca (mh 1002) = (03) 
Vo(Gu ) AIAR) dla fazy nematycznej 


gdzie gr jest cHS w przybliżeniu Percusa-Yevicka /133,63/, 
V jest objętościę wykluczona w przybliżeniu Berne-Pechuckasa 
/132/,a 7 pewną efektywną gęstością. Mulder i Frenkel /126/ 
obliczyli F* za pomoca y-expansion /117/. Sin Doo Lee /130/ 
zastosował wzór (6.1) z g( 9 vd = A4(9)/4% 
Powstały również prace poświęcone badaniu przejścia fazowego 
N-SmA w układzie idealnie uporządkowanych sferocylindrów 
/134-136/. Mydele Muldera /134/ oraz lena i Meyera /135/ zostały 
oparte na rozwinięciu wirialnym gg, ModelSómozy i Tarazony 
/136/ stanowi dosyć skomplikowane uogólnienie modelu Sin-Doo Lee 
na przypadek idealnie uporządkowanych cząstek, 

Teorię SDA zaproponowaną w rozdziale 4 zastosujemy do zba- 


dania przejścia fazowego N-I w układzie twardych sferocylindrów 
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i elipsoid oraz przejścia fazowego N-SmA w układzie twardych 
sferocylindrów o pełnej swobodzie rotacji. To ostatnie przejś- 
cie fazowe nie zostało do tej pory zbadane i opisane w żadnym 
modelu twardych cząstek. Zbadamy także idealnie uporządkowany 
orientacyjnie układ twardych sferocylindrów /brak swobody 
rotacji/. 

Badanie przejść fazowych w ciekłych kryształach jest 
testem dla naszej teorii. Aby poprawnie opisywać zjawiska po- 
wierzchniowe, teoria musi dobrze opisywać zjawiska objętościowe 
a przede wszystkim tworzenie się faz ciekłokrystalicznych. 


W rozdziale przedstawimy wyniki prac /137-139/. 


6.2. Przejście fazowe faza izotropowa-faza nematyczna 
Badanie przejścia fazowego N-I rozpoczniemy od zbadania 
stabilności fazy izotropowej względem zaburzeń o symetrii nema- 


tycznej 
Sgtne)= 93, Au Å (m8), (6.4) 


W modelu Onsagera zagadnienie to było badane /18,140-141/. 
Funkcja korelacji fazy izotropowej obliczona według wzoru 


(4.37) jest równa 


KTG? = 2 Prag (Ø)W (Gefu) + Ag (9):9 [dys dada 

*W (5-13 1684, 093) W(83-%10040) lam)", (65) 
gdzie AY R jest dana wzorami (4.301) i (4.35). Zastosujmy 
wzór (2.39) określający stabilność fazy ciekłokrystalicznej. 
Po prostych przekształceniach wykorzystujących ortogonalność 
wielomianów Legendre'a otrzymujemy warunek stabilności fazy 


izotropowej względem zaburzeń nematycznych. 


http://rcin.org.pl 


er 
n+ Takele) lp (dø, laing R (08,)>0 (6.6) 
= p Ja aid R 2) 


dla 1=2,4:.4+ 

farunek ten jest bardzo podobny do warunku stabilności w modelu 
Onsagera /18,20/. Różnica występuje jedynie w czynniku gęstoś- 
ciowym > „ który w granicy Onsagera jest równy jedności, 
Warunek (6. 5) przestaje być spełniony najpierw dla 1=2. 


Otrzymujemy wówczas warunek na bifurkację N-I w naszyn modelu 


1— TE Deae = 0, (6.7) 
gii 4f 8 (T 4 T) opr. (6.8) 


Równanie (5.7) rozwiązujemy numerycznie metodą Newtona. W ta- 

E 
beli (6.1) przedstawiona zostałagęstość bifurkacji sz = Gaa 
w funkcji L/D. Jest to malejąca funkcja L/D i w granicy 


Onsagera /§4.2/ dęży do zera. 


ı L/D i 10 ı 8 1 6 3 5 14 4 ! 3 ! 2.46 I 
mAr | ag 0.332 | 0.374! _0.4! 0.42810.466 I 0.492 | 
tp ! OJ 1 | 0.51 0.3! 0.28! 0.2 | 
K 

NI | 0.52 ! 0.559! 0.647! 0.714! 0.77910.787! 0.823 I 


Tabela (6.1) 


Dla L/D=0.28 mår przekracza gęstość maksymalnego upakowania 
sferocylindrów N cp” TT (3L/D+2) /(/ŚL/D+YZ)6, co oznacza 
że dla L/D £ 0.28 bifurkacja zachodzi w obszarze niefizycznyn 
Istnieje więc minimalna anizotropia kształtu cząstki, przy 
której może utworzyć się faza nematyczna. W szczególności 


twarde kule nie tworzę fazy nematycznej. 
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W układzie jednorcænym energia swobodna F wyraża się 
wzorem /54.4/ 


AET- fdo St lni) H+ AC), (68) 
dodw Wolwe) (6.10) 
jdodw'Vo(ww ) ir ir 

Jeśli przybliżymy f(w) =f(0) funkcją próbną 

F(0)= exp (LF (100) /anjdownGexp (d Bos), (6.11) 
gdzie dl jest parametrem tak dobranym aby energia swobodna 


F byta minimalna tj. 


6.12 
27 =0 ) (BJ) 


to otrzymamy równanie wiążące 1 id. 


LIŚR (m6) -Q] +20) KROWY =0, (613) 
KR (c056,)-QX)=|dayda, Voleu wa) f(O) (81) > 


«(BC$)-Q) / fadon olco Jin), (6.14) 
T 
< PY(00)) = 2T [dosing £10) R0) , (6.45) 
(o) 


Q jest nematycznym parametrem porządku (1.1), 

Równanie to ma dwa rozwiązania: izotropowe ol =0 i anizotropo- 
we & 40. Aby określić współistnienie N-I uzupełniamy równania 
(86.13-56.15) warunkami na równość ciśnień i potencjałów che- 
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micznych obliczamy ze wzorów 


ir "8 29-13 LA) ) (616) 
A= $ + F. (6.47) 


W ten sposób otrzymujemy trzy równania na trzy niewiadome: 
I» NN 3 L . Tabela (6.2) zawiera otrzymane wyniki tzn, 


ciśńienie współistniejących faz, paranetr porządku fazy nema- 


tycznej we e almen) msme ci NI > NN obu faz oraz skok 


ka JF 


gęstości Am= 


3 ' 2. 
14 ! 2.784! 0.575 !0.407 10.437 | 0.0706, 3.509 | 0.8787 
5 I 2.875! 0.588 10.376 !0.408 ! 0.0809! 2.889! 0.8835 
! 6 ! 2.971! 0.601 10.350 ;0.382 | 0.0897] 2.477, 0.8870! 
7 ! 3.059! 0.614 10.327 10.360 ! 0.0973! 2.176! 0.8898 
! 8 ! 3.168! 0.627 10.306 10.339 | 0.104 | 1.943 į 0.891 
r 9 I 3.2651 0.639 10.287 10.320 !0.110 ! 1.755! 0.8931 
ilo ! 3.361! 0.550 !0.270 10.303 | 0.116 | 1.599 į 0.8943 
f100 | 5.114! 0.783 10.033410.04161 0.219 i 0.1491 0.9056! 
100 t 5.481! 0.799 ! 0 ! 0 10.243 | O | 0.9069! 
Tabela 6.2 


Porównując wyniki umieszczone w tabeli (6.1) i (6.2) zauważamy, 

że linia punktów bifurkacji leży pomiędzy linią 4 n (+2) 
Nil). W granicy L/D->00 odtwarzamy wyniki modelu 

Onsagera. /patrz(4.21/. Niewielkie różnice / ~- 7%/ pomiędzy 


wynikami naszego modelu a (4.21) wynikają z różnej postaci 
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f (w) /por.6.11) i 4.20) użytej w obliczeniach. Punkt L/D=2.46 
jest szczególny, VI $ 6.3.2 pokażemy, że linia bifurkacji N-SmA 
łęczy się w tym punkcie z linią N „(LD . Oznacza to, że dla 
L/D X 2.46 faza nematyczna jest metastabilna a globalne minimum 
F odpowiada w zależności od gęstości bądź fazie izotropowej 
bądź fazie snektycznej. 

Mechanizm przejścia fazowego N-I w naszym modelu jest następu- 
jący: cząstki uporządkowane orientacyjnie "słabiej" ze sobą 
oddziałują niż cząstki nieuporządkowane. Efektywna gęstość M 
jest: mniejsza niż gęstość rzeczywista n i stąd Ate (T) Akal) 
co oznacza że uporządkowanie orientacyjne zmniejsza F, Z dru- 
giej strony część idealna Fi’ /rotacyjna/ zawsze preferuje 
fazę izotropową. Przy pewnej gęstości F zaczyna dominować 


i następuje przejście fazowe N-I. 


6.3. Przejście fazowe faza nomatyczna - faza smektyczna A 
6.3.1. Przypadek idealnego uporządkowania orientacyjnego 


Rozważmy układ równoległych sferocylindrów. W tyn przy- 
padku f ref jest dane wzorem (4.30ii), ponieważ w nieskończenie 
silnym polu zewnętrznym Vext ( W) /S 5.5/, które porządkuje częst- 
ki wzdłuż określonego kierunku w przestrzeni, faza izotropowa 
nie jest stabilna. Rolę układu odniesienia względem którego 
oblicza się energię swobodną stanu niejednorodnego przejmuje 
idealnie uporządkowana faza nematyczna. 

Zbadajmy utratę stabilności fazy nematycznej względem 
zaburzeń o symetrii smektycznej 


DO 


S9(50)=2_ gn toś(nkz) (6.18) 


n=4 
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TA 


które opisują periodyczną zmianę gęstości w kierunku osi zl|n. 
9 n 58 Paranetrami porządku a d=2][/k - okresem smektycznym, 
Ze wzoru (4.37) otrzymujemy gan dla idealnie uporządkowanej 


fazy nematycznej 
— „dN ll 
-kgr c3 (5-5)= 2AVe (F)W%- 1) + Afryle)9* 
x [dr Wa -13)W (1-1). (6.19) 
Waga w (4.33) nie zależy od orientacji w przypadku idealnego 


uporządkowania. Warunek stabilności (2.39) zapisuje w nastę- 


pujęcej postaci 
på lonl {E +24ptø)Mnk) + PA (gó) 0, (6.20) 


gdzie 


W(nk)= faronkz W(T) (6.24) 
jest transformatę Fouriera wagi w. Po prostych obliczeniach 
otrzymujemy 
X (mk) =| pp {sim (nk(D+L)) ~ sin (nkL)} (6.22) 


- MD, wos(nk(D+L))]/(2MD*L+4MV3) ; 


w granicy L/D > eo „ KL=const 


W(nk)= Aln.(nkL) (6.23) 


Bezpośredni rachunek pokazuje, że wzór (6.23) otrzymujemy 


także dla układu idealnie uporządkowanych cylindrów o długości 


L. Oznacza to, że w tym przypadku układ idealnie uporządkowanych 


sferocylindrów o długości [/D>>%0 jest równoważny układowi 
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9. 
cylindröw o dowolnej dłupści Le 


Warunek (6.20) załamuje się najpierw dla n=1 ti. g n=0 n>2 i 


A+ 20069 Abe (e) + K AYA (p)=0. (6.24) 


Gęstość bifurkacji på odpowiada minimum funkcji P(x) spełnia- 
jęcej równanie .(6.24) . Różniczkując (6.24) względem k otrzy- 


mujemy że k* odpowiadające minimum (k) spełnia równanie 
r 
K) 
a. str (6.25) 


Tabela (6.3) przedstawia otrzynane wyniki dla L/D € (0,00) é 
Gęstość bifurkacji jest przedstawiona zarówno w jednostkach 


; 4 
9 cp jak i objętości cząstki v,. 


Tabela 6.3 


Gęstość bifurkacji maleje wraz ze wzrostem L/D, szczególnie 
szybko dla L/D € (0,2). Dla L/D> 5 gęstość bifurkacji jest 
praktycznie stała. Okres smektyczny liczony w jednostkach L 
również maleje gdy L/D rośnie. W praktyce d jest trochę większe 
niż długość cząsteczki L+D. Gdy L/D>0 ytedy d/L> D/L. 


Gęstość fazy smektycznej wyraża się wzorem 
p(2)= 99(kz)=9[1+ 2 Oncostnkz)) (6.26) 


Wstawiając (6.26) do (4.27)-4.29) otrzymujemy energię swobodną 
na jedną cząsteczkę. Nieskończoności wynikające z Vext(60)-1 
Inf (w) wzajemnie się znoszą tak jak to pokazaliśmy w § 5.5 


przy omawianiu idealnego uporządkowania w pobliżu ściany. 
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Korzystając z tego faktu „las F; 


F/Nkgl = Pa 1) +4 d5 l$) inp) t 
+ dr fax 20) (FO) | (6.22) 
9(5)=4 +S on W(nk) cos(mE) (6.24) 


Zminimalizujmy F względem Va i k. Po prostym rachunku otrzy- 


mujemy 


na ję (a crt ye) + Ar J43 mnd Ay 


saso OR) Alons = =0, (6.29) 
dni) _ 


23/Nkgl 1 cosh 
4 fodzo(e) Abe (NPR): >» ika 
SOK > 7 hlo PEAP, Ej 


Równania te mają dwa rozwiązania: nematyczne p „.kjo oraz 
smektyczne nK40. Przeprowadźmy analizę rozwiązań równań 
(6.29) , (6.30) v pobliżu punktu utraty stabilności fazy nematycz- 
nej /analiza bifurkacji/. Pozwoli to ustalić rodzaj przejścia 
N-SmA. W tym celu sprametryzujmy rozwiązanie równań (5.29) (6.30) 


w pobliżu (RZA w następujący sposób 


RY | 
=h tier PEt., 
= k* + Ke + BĘ. (6.30) 
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Wstawmy (6.31) do (6.29) i (6.30) . Ponieważ równania te muszą 


być MAPY dla dowolnego € więc 


dm KU 
den! (RE 1 


Dla m=12odtwarzamy warunki bifurkacji (6.24). Dalsze rachunki 


(6.32) 
m=4,2.... 


prowadzą do następujących wyników 
a> Øre 
$2= EE + OlE?), 
Pn=O(e") , 
ngt PE ole) ; 
k=k* + OE). 


Zaznaczyliśmy w jawnej formie tylko te współczynniki rozwinięć, 


(6.33) 


które są niezbędne do określenia rodzaju przejścia fazowego. 

2 2 
Współczynniki sj i 7” są podane w dodatku C. ør wyznacza 
skalę. Eliminujac & z równania na 8 1 wstawiając je do 


równania na 94 otrzymujeny 


=p" (6.34) 


Zostało sprawdzone numerycznie, że 4, > O a więc rozwiązanie 
smektyczne istnieje powyżej punktu bifurkacji (1>7) . Oznacza 
to, że przejście fazowe jest ciągłe i że punkt bifurkacji n* 
jest punktem przejścia fazowego, Q2 A P> są również rosnącymi 
funkcjami N tak samo jak dominujący parametr porządku 94 . 
Zanik Pi przy nd jest dany przez zależność potęgową z wykład 


nikiem 1/2 charakterystycznym dla teorii pola średniego /142/. 
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W pobliżu punktu bifurkacji różnica energii swobodnych fazy 


ae iy fazy nematycznej wyraża się wzorem 


2) cH 

> Jag FI (ge) Ki e + OLE) (6.35) 
Ponieważ żę < 0 i > O więc Fx 7" co potwierdza 
wniosek o rodzaju przejścia fazowego. Rozwięzanie smektyczne 
powyżej punktu n* ma niższą energię swobodną niż rozwiązanie 
nematyczne, aw NĄ = rË obie energie są równe. Punkt N = a 
jest więc punktem przejścia fazowego ciągłego. 

Mechanizm przejścia fazowego N-SmA w układzie idealnie 
uporządkowanych sferocylindrów jest następujący: 
Cząsteczka z fazy smektycznej ma większą swobodę ruchu transla- 
cyjnego w płaszczyźnie xy /płaszczyzna równoległa do warstw 
smektycznych/ i oddziałuje z mniejszą liczbą cząstek niż cząs- 
teczka z fazy nematycznej. Zysk entropowy związany zutworzeniem 
warstw opisuje trzeci wyraz we wzorze (6.27) wzięty ze znakiem 
przeciwnym. Z drugiej strony ruch danej cząstki w fazie smek- 
tycznej jest ograniczony do jednej warstwy. Zmniejszenie entro- 
pii translacyjnej związane z tym efektem jest opisywane przez 
drugi wyraz we wzorze (6.27) wzięty ze znakiem ujemnym. 
Przy pewnej gęstości wkład wyrazu trzeciego we wzorze (6.27) 
do energii swobodnej staje się dominujący i wtedy następuje 
przejście fazowe. Zrozumiałe jest, że gęstość przejścia fazo- 
wego n* maleje z L/D. Im dłuższe są cząstki tym większy 


jest zysk entropowy związany z utworzeniem warstw. 
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6.3.2. Przypadek nieidealnego uporządkowania orientacyjnego 
Zbadajmy teraz układ sferocylindrów w fazie nematycznej 
bez pola zewnętrznego. Funkcja korelacji cz /patrz (4.377 


wyraża się wzorem 


skra (6-5.,0060,)=2A%g (FIN (12100) + (5-1) * 
Aeg (P) AP) Alu) + A Vg (P) AP fdo Fflos) 
xW (22160403) + ABe (5) A Pl) fde $ (03) har 08) 
+ Aeg (9) gf dodwyd WT 04105) W (5-5, 04): 
* 9*4(03) lwy), (6.36) 
gdzie P jest dane wzoren (6.12) a f (©) wzorem (6.11). cz” ma 


symetrię fazy nematycznej. 


Zbadajmy utratę stabilności fazy nenatycznej względem zaburzeń 
o symetrii smektycznej re 
z= nK 
Setru)= flo)? Pn cosnk (O he S~ (6.37) 
h=4 n= 
Stosując wzór (2.39) i korzystając z wyrażenia 
f do, $ Cor) AP) =D, (6.38) 


otrzymujemy warunek na stabilność fazy nematycznej względem 


zaburzeń snektycznych, 


på (al? f 4+ 29 Veg (P) wy (nk sd) 4 dr de 


k 
Wy(nkid) = |w (X. Oy JE tru) fl) drduydn, 
http://rcin.org.pl (6. 40) 
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gdzie o. jest dane równaniem (6.13). Zakładamy, że warunek 
(6.39) załamuje się najpierw dla n=3 


A + 29 Ag (P) Wu (Kid) + Abe (PWy(Kia)=0. (6.14) 


Wektor k dla którego (5.41) jest spełnione dla minimalnej 


gęstości P jest rozwiązanien równania 
das Kd) o, (6.42) 


Ukžad równań (6.42) (6.41) i (6.13) został rozwiązany numerycz- 
nie. Cażki siedmiokrotne występujące w (6.40) zostały zreduko- 
wane do całeK pięciokrotnych przez całkowanie względem c z wy- 
korzystaniem parametryzacji bryły wykluczonej objętości podanej 
w 8 4.1. 

Tabela (6.4) zawiera uzyskane wyniki tj. gęstość bifurkacji NN ap 


I L/D | KL d* vat PV, d/L | dAL+D) | a ! 


3  |14.499| 3.668 | 0.496 14.867 | 1.397 | 1.047 10.682 


t i 2i ' 1 i 1 t i 
————4————mm ZZ ZZ ZZ 
! 4.4961 8.477 1 0.528 14.042 ! 1.398 ! 1.165 10.876 


5 i 
i 
6 1 4.525 111.829 |, 0.553 14.034  ; 1.389 į 1.190 ;0.913 H 
8 i 
1 
l 


I 10 14.622 130.186 | 0.641 14.310 11.359 | 1.236 


Tabela (6.4) 
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="B1 = 
wektor Falówy og paranetr d* „ parametr porządku nematycznego 
Q, ciśnienie p oraz okres smektyczny d. Porównując tabelę (6.4) 
z tabelą (6.2) zauważamy, że dla L/D=2.46 linia bifurkacji 
N-SmA łączy się z linią Ny (1/0). Pznacza to, że dla L/DX 2.46 
następuje bezpośrednie przejście fazowe faza izotropowa =- faza 
smektyczna. Gęstość bifurkacji szybko rośnie z L/D dla L/D 23, 
maleje zaś dla L/D £ 3. Parametr porządku rośnie wzdłuż 
całej linii bifurkacji. Wyniki te szczegółowo przedyskutujemy 
po przeprowadzeniu analizy bifurkacji. 

Przybliżamy lzw) dla fazy smektycznej następującym 


wzorem: 


glo) =P lu) (A+ Å fncosnke )= pflwjp(kz). (6.43) 


Przybliżenie to spełnia warunek periodyczności (2.12) niemniej 
pomijamy w nim zależność orientacyjnej funkcji rozkładu od 
zmiennej z. Energia swobodna na jedną cząsteczkę jest równa 


/patrz (4.297 


217 
F/NkgT = ż,|do(dz Ho) pE) lm LOA), -4f+ 
2r 
+ dr Jd3 PE) A ht (175), (6.44) 
PE)=Wy (0,4) +Ż gn Wy (nkid)cosn$ „ (6.45) 


Dla P nfo jest to energia swobodna fazy smektycznej, zaś dła ' 

Q n=0 (6.44), redukuje się do wzoru (6.9) na energię swobodną o 
fazy nematycznej. 

Minimalizujeny Ff ze względu na Pr’ Lik. w wyniku otrzymuje- 


my następujące równania: 
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Get Å fa d% osn% In PE) + frag coon Ahe (NPE) 

+ år faz ę6)t, 1 (PEJN kA) coan = =0, (646) 

3009. 4 år falsk TON pen nA) 0 
(eg 


n=4 


OSA | (CBD -Q) dg y(5 AW, (PJ) 
rn KLACZ 2 MCO pin ) =0. (6.48) 


Parametryzację rozwiązań równań (6.46) - (6.48) w okolicach 
punktu bifurkacji przeprowadzamy tak jak w § 6.3.1 i badamy 
kolejne rzędy rozwinięcia równań w € .w zerowym rzędzie w € 
odtwarzamy równanie (6.13) „aw pierwszym rzędzie w & rów- 
nanie bifurkacji (6.41) . W drugim rzędzie w £odtwar :any równanie 


(6.42) oraz otrzymujemy 


L= q"=O, i (6-49) 
POV BE RAM re. å gå (650 


GEN | td gy 


Trzeci rząd Syrię w E prowadzi m.in. do równania 


2 ) 90 PG 3 2'F 
EP er 2) eo (6.52) 
2 Å [2 OT pl 


alogi JA o +2grf" ole 0 


2 Fl = 
areta (650 
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Z równania (6.51) i (6.52) obliczamy dim, z dokładnością 
do z” różnica energii swobodnej fazy smektycznej i nema= 


tycznej jest równa 


1) jg! 
AF" PPE (63) 


PE. Er (2 Nn am z równań otrzymujemy 


peia T (6.51) 


PLAZE + ŻĘ” zg Grim) (6.55) 


Z obliczeń numerycznych wynika, że 


De | + ŻE. (9 56 
Zlec” * Bnk (6.56) 


więc znak wf) decyduje o rodzaju przejścia fazowego. Dla 
L/D > 3.29 1/2 O i przejście jest ciągłe, zaś dla L/D< 3.29 
1/7 L O i przejście jest pierwszego rodzaju. W tym drugim 


AF DO, qarisa (6.57) 


co oznacza, że istnieje "gałąź" rozwiązań smektycznych o więk- 


szej energii niż energia fazy nematycznej w pobliżu LOT ZĘ 
= * 


Vystępowanie takiej "gałęzi" rozwiązań jest charakterystyczne 
dła przejść fazowych pierwszego rodzaju /141/. Dla L/D= 3.29 
linia przejśćkwpierwszego rodzaju łęczy się z linią przejścia 
ciągłego. Punkt ten nazywamy tröjkrytycznym. 

Mechanizm tworzenia fazy smektycznej w omawianym układzie 
jest związany,tak jak w przypadku układu cząstek równoległych 


ze zwiększeniem entropii translacyjnej w płaszczyźnie warstw. 
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Zysk ten jest większy niż zmniejszenie entropii translacyjnej 
związanej z ruchem cząstek w kierunku prostopadłym do warstw, 
Taki scenariusz nie jest jednak kompletny, ponieważ w układzie 
z orientującym polem zewnętrznym /idealne uporządkowanie/ 
MN -sma 7/0) maleje z L/D, a w układzie bez pola wielkość 

ta rośnie z L/D. Aby wytłumaczyć te różnice zbadajmy wielkość 
x=L {sin wzdłuż linii przejścia fazowego N-SmA. Wielkość 
ta opisuje średnie odchylenie cząstki w kierunku poprzecznym 
do osi z. Aby utworzenie warstw prowadziło do zwiększenia 
entropii translacyjnej cząstki nuszą być silnie uporzędkowane 


a więc x musi być małe. Tabela (6.5) pokazuje x w funkcji L/D. 


! L/D 
1 g/h 
1 


i 
x/D | 


4 24.5 |5 16 48 140 
1 
I 


' 
1 
1.256! 1.259! 1.260! 1.267 1.309i 1.324 


Tabela (6.5) 


Wielkość ta bardzo wolno rośnie wzdłuż linii przejścia fazowego 

w przeciwieństwie do parametru porządku Q / tabela (5.4)/. W przy- 
padku idealnego uporządkowania pole zewnętrzne wymusza x=0. 

W przypadku bez pola zmniejszenie x można wymusić tylko zwiększa- 
jęc gęstość. Dlatego EN (1/0) rośnie z L/D i oczywiście 
jest większe niz w przypadku idealnego uporządkowania. 


6.4. Utrata stabilności fazy izotropowej względem zaburzeń 
smektycznyc 


Zaburzenie względem, którego badamy łamanie stabilności 


fazy izotropowej ma teraz następującą postać: 


ó9(010) = A a nu £h(nkr) R (w) (6.58) 
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Faza smektyczna jest uporządkowana zarówno orientacyjnie jak i 
translacyjnie, stąd mieszany charakter zaburzenia (6.58 ) /19/. 
Wzory (6.5) na og oraz (2.39) prowadzą do warunku stabilności 


fazy izotropowej względem zaburzeń (6.58). 


5 © uga) de Bla) Pa Coon) + 2 AV)» 


(6.59) 
x [drday dw, WOT oy, z) OS (MK) Ri (0) Å (car) + 


+ Are (9) 9 -fda deog dip, WCG ca, 094 )eøs (nk): 


* Ja den FR P (0) (110005) wslnkrz)?>0 


Dla każdego n (6.59) musi być spełnione tzn. 

2 Nu One Qh >0 (6.60) 
gdzie W,, , jest dane przez wyrażenie w nawiasie klamrowym z 
(6.59) . 
Zakładamy, że stabilność załamuje się najpierw dla L,L'=0,2. 
W tym przypadku aby zapewnić dodatniość formy (6.60) muszą być 
spełnione trzy warunki /19/. 


Woo >O, (6.64) 
W22>0, (6.62) 
WooWzz — Woo >O (6.63) 


Zerowanie się któregokolwiek z wyrażeń (5.61) - (6.63)wyznacza 
pewną linię AU). Linia utraty stabilności jest wyznaczona 
przez te gęstości w; które są najmniejsze dla danego L/D. 


Obliczenia numeryczne zostały przeprowadzone w przedziale 
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(0,2.46) ponieważ w tym przedziale występuje bezpośrednie 
przejście fazowe faza izotropowa-faza smektyczna A /patrz 3 6.3/. 
Dla 2 < L/D< 2.46 najpierw warunek (6.62) przechodzi w rów- 
ność przy czym bifurkacja następuje dla k=0. Oznacza to, że 
pomimo, iż spodziewamy się przejścia fazowego I-SmA bifurkacja 
zachodzi dla zaburzeń czysto nematycznych. W przedziale 1.5< 
L/D< 2 linię bifurkacji wyznacza warunek (6-63) /kfo/. Dla L/DC 15 
stabilność łamie się dla zaburzeń czyśto smektycznych, a linię 
bifurkacji wyznacza zerowanie się 6.61). Wyniki liczbowe 


zostaną przedstawione na diagramie fazowym w 8 6.6. 


6.5. Punkt potrójny: Współistnienie SmB-SmA-N 


Smektyk B jest scharakteryzowany periodycznyn bądź quasi- 
periodycznym porządkiem heksagonalnym w warstwach. Załóżmy, że 
krystalizacja w warstwach w przejściu fazowym SmA-SmB zachodzi 
podobnie do krystalizacji dwuwymiarowego układu twardych dysków. 
Jeżeli tak jest to przejście fazowe SmA-SmB zachodzi przy takiej 
gęstości cząstek w warstwie przy której krystalizują twarde 
dyski. Proste przeliczenie gęstości powierzchniowej Ps krysta- 
lizacji dysków na gęstość objętościową prowadzi do wzoru 


na gęstość przy której tworzy się SmB, 


p= 05 re (6.64) 
gdzie d jest okresem smektycznym. Wielkość d jest ściśle 
policzona na linii przejścia fazowego N-SmA. Twarde dyski 
krystalizują przy gęstości 4 9 s/ T D?=0.69 /143/. Obliczenie 
wykorzystujące wzór (6.64) pokazuje źe dla L/D=6 fN-GnA" Yo" 
=0.55 /d=1.3895/. Punkt ten możemy uznać za punkt potrójny, 
albowiem współistnieją w nim trzy fazy: nematyczna, smektyczna A 


i smektyczna B. 
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6.6. Diagramy fazowe, porównanie z symulacjami komputerowymi 
i innymi teoriami 


6.6.1. Porównanie dla układu sferocylindrów równoległych 


Rysunki (6.1.) i (6.2) OE JED odpowiednio gęstość 
przy której zachodzi przejście fazowe N-SmA i okres smektyczny 
w punkcie przejścia jako funkcje L/D. Linią przerywaną zostały 
zaznaczone wyniki symulacji. W symulacjach komputerowych er 
zmienia się od 0.39 dla L/Deo do 0.56 dla L/D=0.5. Nasze 
wyniki /tabela (6.3) i wykres 6.17 sa zgodne z symulacjami 
w granicach 25%. Okres smektyczny rys. (6.2)zgadza się z wynikami 


symulacji z dokładnością 10-cio procentową. 


rys 6.1) gęstość pre pain fazowego N-SmA w funkcji długości 


sferocylindra dla układu równoległych sferocylindröw.Linia 
przerywaną został zaznaczony wynik symulacji komputerowych. 
Linią ciągłą - wynik modelu 
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rys .6.2) okres snakiycziy w funkcji długości sferocylindra. 


Linia przerywana - wyni 


symulacji. Linia ciągła - wynik 
modelu. 
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W symulacjach przejście fazowe N-SmA obserwuje się aż do 
L/D=0.25. Dla L/D < 0.25 następuje bezpośrednie przejście fazo- 
we ciecz-ciało stałe. W naszym modelu faza smektyczna występuje 
aż do L/D=0 tj. nawet dla twardych kul. Dla L/D=0 zostało zka- 
dane przejście fazowe ciecz-ciało stałe w tym modelu /51/. 
Pin (L/D=0) jest mniejsza niż 
gęstość płynu twardych kul współistniejących z ciałem stałym. 


Okazuje się, że gęstość 


Ten niefizyczny wynik jest spowodowany uproszczoną postacią wagi, 
która w przybliżeniu zaproponowanym w rozdziale 4. nie zależy 
od gęstości, L/D=0.25 wyznacza granicę stosowalności zapropono- 
wanego modelu w układzie idealnie uporządkowanych sferocylindröw. 
Wyniki zgodne z symulacjami z dokładnością do 10% w całym 
zakresie L/D A Somoza i Tarazona /136/. Warto podkreślić, 
że krzywa g (1/0) w ich pracy ma identyczny kształt * ta 
uzyskana w modelu SDA przedstawionym w rozdziale ar AA 
tylko przesunięta jako całość w stronę wyższych gęstości, 
Dzięki temu, że g (10) szybko rośnie gdy L/D- O model Somozy 
i Tarazony przewiduje przejście fazowe ciecz-ciało stałe przed 
pojawieniem się fazy smektycznej. Mulder /134/ zbadał bifurkację 
w swoim modelu w granicznym przypadku L/D=oo i uzyskał wynik 
GU $p=0.41 zgodny z synullcjami. Autor jednak nie pokazał, że 
przejście fazowe N-SmA jest ciągłe, Z kolei w modelu ena i 
Meyera /135/ uzyskuue się ciągłe przejście N-SmA ale dla gęs- 
tości de =0.2 niezależnej od L/D. 
Podsumowując: jakościowo jak również ilościowo zapropono- 
vany model dobrze opisuje zjawisko przejścia fazowego N-SmA 


w układzie równoległych sferocylindrów o długości L/D > 0.25. 
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6.6.2 Porównanie dla układu sferocylindrów nierównoległych 


Rysunek (6.3) przedstawia globalny diagram fazowy układu 


twardych sferocylindrów bez orientującego pola zewnętrznego. 


0.7 


0.2 


rys.(6.3) Diagram fazowy układu sferocylindrów 

I - faza izotropowa, N- faza nematyczna, SmA ~ faza smektyczna A 
A- wyniki symulacji, - punkt pA ai Å - punkt potrójny 
/N-SmA-SmB/, linia przerywana - linia bifurkacji, linia ciągła = 
linia przejścia fazowego ; 
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Zawiera on fazę izotropową, I, nematyczna,N, i smektycznę A,SmA. 
Linie przerywane oznaczają linie bifurkacji a linie ciągłe to 
linie przejść fazowych. Zakropkowany obszar odpowiada duufazowemu 
obszarowi N-I. Czarny punkt oznacza punkt trójkrytyczny a krzyż 
to hipotetyczny punkt potrójny tj. punkt współistnienia N-SmA-SmB, 
Białe trójkąty są wynikami symulacji komputerowych. U symulacjach 
zmierzono tylko 9 1Vo " przejściu fazowym N-I a nie zmierzono 

9 Nio; Stad tylko 1 punkt zaznaczony na wykresie dla przejścia 
fazowego N-I. Symulacje komputerowe zostały przeprowadzone dla 
jednej długóści sferocylindra L/D=5./415/ ` 

Punkt trójkrytyczny leży bardzo blisko linii przejścia N-I. 
Jest to charakterystyczne dla rzeczywistych układów ciekłokrys- 
talicznych. Zmiana przejścia fazowego z cięgłego na pierwszego 
rodzaju spowodowana jest silnym sprzężeniem pomiędzy nenatycz- 
nym parametrem porządku Q a smektycznym ọ 1° Dla dużych L/D 
Q jest bliskie jedności i sprzężenie z Pi jest zaniedbywalne. 
Wówczas przejście fazowe jest ciągłe. Gdy linia przejścia M-SmA 
zbliża się do linii przejścia N-I Q silnie sprzęga się z Pa i 
przejście staje się pierwszego rodzaju /2/. 

Z diagramu (6.3) możemy odczytać praktyczny zakres stosowal- 
ności modelu SDA /rozdział 4/ w układzie sferocylindrów nierów- 
noległych.- Dla L/D=0 otrzymujemy taki sam niefizyczny wynik 
jak w przypadku idealnego uporządkowania /sferocylindry równo- 
ległe/. Zauważany również, że gęstość bifurkacji 1-SmA szybko 
maleje gdy L/D zmienia się od 1.5 do O. Z fizycznego punktu 
widzenia gęstość bifurkacji powinna raczej rosnąć,gdy L/D 
malejeytak by mogła przeciąć linię przejścia fazowego faza 


izotropowa - ciało stałe dla małych L/D. 
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Sędzę, że SDA poprawnie opisuje układ sferocylindrów dla 

L/D > 1.5. 

Vykres (6.4) przedstawia nematyczny parametr porządku w funkcji 
L/D wzdłuż linii przejść fazowych, a rysunek (6.5) okres snek- 
tyczny w funkcji L/D. Zmniejszanie się okresu smektycznego 

wraz z L/D związane jest ze zmianą rodzaju przejścia fazowego. 
Dla dużych L/D okres smektyczny zachowuje się tak jak w przy- 


padku idealnego uporządkowania tzn. wolno maleje z rosnącym L/D. 


1.0 


Q 


0.75 


rys 6.4) nenatyczny parametr porządku w funkcji 
długości sferocylindra wzdłuż linii przejść fazowych 


/patrz rys (6.3 
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1.45 


135 


rys {6.5) okres smektyczny w funkcji długości sferocylindra 
wzdłuż linii przejścią fazowego N-SnA /patrz rys. ©.3y 


W tabeli (6.6) przedstawione jest porównanie wyników SDA oraz 


wyników pracy Sin Doo Lee /144/ dla przejścia fazowego N-I, 


L/D=5. 
pp O e O -E 
i L/D=5 pa Po 1279 Q 
i i 1 1 kn i ' 
isymulacje i i H i 
Ikomputerowe | 0.4 JRE ! 4,9 10.3-0.4 
! 125 i 1 i [| i 
ISin Doo Lee i ær t 1 1 I 
1/144/ | 0.399 1 0.417 15.36 10.67 | 
TESTENE TEJ GT" EET TW EEE 
SDA /438/ i 0.38 I 0.41 12.9 10.59 I 
s A) 


Tabela (6.6) 
Wyniki obu teorii zgadzają się z wynikami symulacji komputerowych. 
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6.6.3 Porównanie dla układu elipsoid obrotowych 

Wiele prac teoretycznych zostało poświęconych przejściu 
fazowemu N-I w układzie twardych elipsoid obrotowych o pół- 
osiach a,b=c/126-130/. Również dla tego układu wykonano wiele 
symulacji komputerowych /119-120/ i dlatego zastosujemy SDA 
do zbadania tego układu. 

Objętość wykluczoną (w Wa) dwóch elipsoid obroto- 
wych można obliczyć rozwijając ję w szereg wielonianów 


Legendre'a /ł45/ tj. 
2L+4 Faux) A XM 
Volonk)edvet 2 Mr Latol (Gr) 84 | ś5 (6.65) 


gdzie x=a/b. Dla x 1 otrzymujemy układ bardzo długich wąs- 
kich cząstek /'"igieł"/ a w odwrotnej granicy x & 1 układ płas- 
kich cząstek. Współczynniki ag (x) i b, (x) oblicza się metodą 
podaną przez Isiharę /145/. Ponieważ /126/ 
(x 4 

ee x ET (6.66) 
więc zaproponowany model SDA przewiduje całkowitą symetrię 
pomiędzy układem elipsoid o stosunku półosi x i 1/x. Wynik 
ten jest zgodny z wynikiem symulacji komputerowych. W symula- 
cjach obserwuje się prawie całkowitą symetrię pomiędzy układem 
elipsoid o stosunku długości półosi x i 1/x. W celu porów- 
nania opisywanej teorii z innymi modelami wykonaliśmy obliczenia 
dla x=3 /1/3/. Współczynniki rozwinięcia zostały obliczone 
/126/ aż do 1=20. Wyniki obliczeń umieściliśmy w tabeli (5.7) 


wraz z wynikami symulacji komputerowych oraz innych modeli. 
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ikomputerove 


Tabela (6.7) 


Uyniki SDA dość dobrze odtwarzają wyniki symulacji komputero= 
wych oraz zgadzaję się z wynikami innych modeli. Należy pod- 
kreślić, że iz wyjątkiem modelu Sin Doo Lee pozostałe modele 
służą do opisu tylko przejścia fazowego N-I w układzie twar- 
dych elipsoid, podczas gdy SDA opisuje znacznie bogatszą 
klasę zjawisk, Wspomniany wcześniej model Sin Doo Lee opisuje 
z kolei tylko układy jednorodne. 

Podsumowujac,teoria SDA dobrze opisuje tworzenie się faz 
ciekłokrystalicznych i może być stosowana w badaniu układów 


niejednorodnych. 
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7. FODSUMOWANIE WYNIKOW PRACY I DYSKUSJA 
7.1 Zestawienie najważniejszych wyników pracy 


1/ W pracy sformułowano teorię SDA /smoothed density 
approximation = przybliżenie wygładzonej gęstości/ dla nie- 
jednorodnych i anizotropowych układów złożonych z twardych 
cząstek. 

2/ Zastosowano SDA do zbadania przejść fazowych w układzie 
twardych sferocylindrów oraz w układzie twardych elipsoid obro- 


towych. 
a. Po raz pierwszy w teorii twardych anizotropowych cząs- 


tek o pełnej swobodzie rotacji zbadano powstawanie fazy smek- 
tycznej A i wyjaśniono mechanizm jej powstawania. Stosując me- 
todę bifurkacji zlokalizowano punkt trójkrytyczny przejścia 
fazowego faza nematyczna =- faza smektyczna A w układzie: twar- 
dych sferocylindrów. 

b, Wyjaśniono również mechanizm powstawania fazy smektycz- 
nej A w układzie sferocylindrów poddanych działaniu nieskończe- 
nie silnego orientującego pola zewnętrznego powodującego idealne 
uporządkowanie. "ykazano, że przejście fazowe faza nematyczna -= 
faza smektyczna A jest w tym układzie drugiego rodzaju. Obliczono 
gęstość przejścia fazowego. 

c. Zbadano bifurkację fazy izotropowej względem zaburzeń 
nematycznych oraz zbadano przejście fazowe faza izotropowa - 
faza nematyczna. Podano nechanizm przejścia fazowego. 

d. Obliczono funkcję korelacji Ornsteina-Zernike dla 
fazy izotropowej i fazy nematycznej. 

e. Zbadano bifurkację fazy izotropowej względem zaburzeń 


smektycznych. 
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f. Wyznaczono diagram fazowy układu sferocylindrów. 

g. Wykazano błędność modelu Tarazony /51/. 

h. Zaproponowano mechanizm tworzenia fazy smektycznej B 
oraz podano wzór na gęstość przejścia fazowego faza smektyczna A 
- faza smektyczna B. Obliczono punkt współistnienia fazy nema- 


tycznej - smektycznej A i smektycznej B /punkt potrójny/. 


3/ Uogólniony model Onsagera /graniczny przypadek SDA/ 
zastosowano do zbadania zjawisk powierzchniowych w ciekłych 
kryształach modelowanych twardymi sferocylindrami i cylindrami. 

ań Wyjaśniono mechanizm porzędkowania orientacyjnego 
cząsteczek ciekłego kryształu przy powierzchni rozdziału fazy 
izotropowej i fazy nematycznej ,N-I, fazy nematycznej i pary,M-V, 
oraz fazy nematycznej i ciała stałego /twarda $ciana/, N-S. 

b. W przybliżeniu ostrej powierzchni zbadano zależność 
napięcia powierzchniowego N-I od kąta pochylenia © t 

c. fl przybliżeniu lokalnym obliczono zależność napięcia 
powierzchniowego N-S, y „ od Og 

d. Obliczono równowagowy kąt or dla powierzchni N-I, 
N-V i N-S. 

e. Pokazano, że napięcie powierzchniowe N-I obliczone ze 
wzoru Kirkwooda-Buffa w przybliżeniu ostrej powierzchni jest 
zgodne ze wzorem obliczonym bezpośrednio z definicji powierz- 
chniowej części wielkiego potencjału termodynamicznego w tym 
samym przybliżeniu. 

f. Obliczono profile parametrów porządku na powierzchni 
N-S w lokalnym przybliżeniu. Zbadano pierwszą nielokalng po~ 
prawkę do przybliżenia lokalnego i przedyskutowano jej wpływ 


na profile, 
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g. Zbadano zależność napięcia powierzchniowego od kształtu 
cząstek w układzie idealnie uporzędkowanym orientacyjnie przy 
twardej ścianie. Pokazano różnice w 4 (04) pomiędzy układem 
cylindrów i sferocylindrów, Podano analityczną zależność y 
od 6, dla cylindrów. 


4/ Sformułowano równanie hydrostatyki dla układów z30z0- 
nych z cząstek anizotropowych w zewnętrznym polu oraz wyprowa- 
dzono ścisłe reguły sum dla ciśnienia, powierzchniowej części 
potencjażu termodynamicznego i adsor5bcji w układzie twardych 


anizotropowych cząstek przy twardej ścianie. 


7.2. Dyskusja 


Przy próbie zastosowania SDA do badania układów ciekło- 
krystalicznych ograniczonych powierzchnią przekonujemy się, że 
równania całkowe na jednocząstkoczkovą zredukowaną funkcję 
rozkładu Qiz) /@.24) patrz też(5.1)/ sa niepraktyczne. Spo- 
wodowane jest to tym, że Pl(zw) jest funkcją aż trzech zmien- 
nych i tym samym jądra całkowe występujące w tych równaniach 
są funkcjami aż sześciu zmiennych. Załóżmy 9(zjw) w postaci 
funkcji próbnej /146/. 

- Peel). gta) At Bosy By 2 ( 1) 
gewo) = © 1 o CET BAJ SP 7. 


gdzie X 1, jest tensorem o symetrii tensora parametru uporząd- 
kowania Q,,. Wzör (7.1) jest uogólnieniem wzoru (6.13) na układy 


niejednorodne, Wstawmy (7.1) do QLA) i zminimalizujmy 
otrzymany funkcjonał względem g(z) i dijlz) e Otrzymujemy: 
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Sy _ Elv S) - 9 
TA = fdzdw EF) E 0 (+ ) 


SA - SQV SE = 


Uproszczenie wyjściowego problemu polega na tym, że teraz 
nieznane funkcje glz) i dC) nie sę już w odröznieniu od 
(zw) funkcjami zmiennych kątowych W = (8,4). 

Równania G.1)- (7.3) stanowią praktyczne uzupełnienie teorii SDA, 

Wiele własności faz ciekłokrystalicznych zostało wyjaśnio- 
nych w oparciu o model twardych, anizotropowych cząstek. Trudno 
jednak odpowiedzieć na pytanie w jakim stopniu efekty wykluczo- 
nej objętości są odpowiedzialne za zjawiska zaciiodzące w ciekłych 
kryształach. Istnieją dwa dodatkowe czynniki: anizotropowo siły 
przyciągające van der ”aalsa oraz giętkie łańcuchy alkilowe 
/122,147/, o których, przy opisie własności ciekłych kryształów, 
należy pamiętać, 

Typowa cząsteczka ciekłego kryształu składa się z dwóch 
/1ub większej liczby/ sztywnych pierścieni benzenowych z "przy- 
czepionyni" łańcuchami alkilowymi /148/. Zdarza się, że jeżeli 
Łańcuch alkilowy w cząsteczce jest zbyt krótki lub nie ma go 
wcale wówczas układ złożony z takich cząsteczek krystalizuje bez- 
pośrednio z fazy izotropowej i nie tworzy faz ciekłokrystalicz- 
nych w odróżnieniu od układu złożonego z cząsteczek o takim sa- 
mym sztywnym rdzeniu ale za to dłuższym łańcuchu alkilowyn. 
Oznacza to, że giętkie łańcuchy stabilizują fazy ciekłokrysta- 
liczne względem ciała stałego; same jednak nie tworzę struktur 


ciekłokrystalicznych. 
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Jak wynika z modelu Maiera-Saupe /95/ w układzie twardych 
kul oddziałujących anizotropowymi siłami przycięgającyni tworzy 
się faza nomatyczna. Z drugiej strony wiadono, że cząsteczki, 
których rdzeń ma symetrię sferyczną oddziałują izotropowyni 
siłami Van der iaalsa, tak więc anizotropia kształtu cząsteczki 
jest niezbędna do wytworzenia anizotropowych sił przyciągających, 

vi ogólności,problem określenia roli różnych czynników 


w tworzeniu struktur ciektokrystalicznych pozostaje otwarty. 
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Dodatek A 


Funkcja AC] 


Przepiszmy wzór (5.5) w następującej,równoważnej postaci: 
A 
Va (core) = fda (Ra) OCR) X (ah 10) (4.4) 


gdzie X jest funkcję charakterystycznę bryły wykluczonej 
objętości /3 4.1/ równą 1 wewnątrz i O na zewnątrz bryły. 

k jest wektorem normalnym do powierzchni a o funkcją 
leaviside'a. W granicy Onsagera /§ 4.2/ możemy pominąć części 
cylindryczne i sferyczne bryły wykluczonej objętości ponieważ 
dają znikomy wkład do całki (A.1) w porównaniu z częścią we- 
wnętrznę O /rys. (4.2) 6.3) /. Skorzystajmy ze wzoru (4.8) 


i zamieńny w (A.1) zmienną [12 na r,s,t. i/ wyniku otrzymujemy 


44 4 
Vi (unon) =+ V jdajdt/arO(rA+0B + ŁC) (rAtoB+ŁC) (AD) 


4 74 A 


gdzie v,=2L?p [sin o 12l jest objętością graniastosłupa; 
A=af, B=bk, C=ck, 
pr) ~ A 
*Iprowadźmy zmienną z=rA+sB+tC, Wtedy 
M 4 AŁSB+ŁC (A3 
V (ac, )= Zb Jas dt | dz OZ. 3) 
TT AAB 
Załóżmy, że AY BXC>0. Mamy teraz do rozważenia dwa przy- 
padki: 1.A> B+C i 2.AZ B+C. 
Przypadek 1. 
Całkowanie w (A.3) rozciąga sie od z=0 do z=A+sB+tC ponieważ 
-A+sB+tC 4 -A+B+C 4 O oraz A+sB+tC >A-B-C20. 
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Vi = Vo (3A%+8*+C*)/42A (4.4) 
Przypadek 2. 


Gdy A4B+C całkę po z obliczamy w następujący sposób: 
AtóB+ŁC 


| dzz0(z)= 4 (A+48+ ŁC) (A+aB+tC) — 
SAW — 1(-A+4B+tCY OleAroBtEC) (A-5) 


Całkujęc pierwszy wyraz prawej strony wzoru (A.5) posit 


otrzymujemy 

aage | (4+8*0)— (A+B-C)- (A-B+C)'] (A6) 
a całkując drugi wyraz dostajemy 

- Am CA+B+C)! (4.2) 


Wstawiając (A.6) i (A.7) do (.3) otrzymujeny Vi w przypadku 2. 
Dla towolnych znaków A,B,C i relacji między nimi oba przypadki 


możemy podsumować w postaci dwóch wzorów: 


_Ay f2[max(IAliBl, ICI)]?+ A%B% C2 
W =4ż% max CIN ETTE j (4-8) 


jeśli lal (Bl. fe] nie spełniają nierówności trójkąta oraz 
Va = Åp Vo CIAL+ Bl+ IE (IA|*1BI-Icl)' - 


(4.3) 
- (IAl-IB| +10)" - Iah 181 +10)" /latølle 


jeśli|A|,|B|,|c] spełniają nierówność trójkąta. 
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Dodatek B 
Obliczanie przekrojów bryły vykluczonej objętości 
Przepiszmy równanie (5.4) w następującej równoważnej postaci 


V(lzzlcywz) = fdr (ga ow.) A Rna- Izl) (8.4) 


Funkcja ta jest równa polu powierzchni przekroju bryły wyklu- 

czonej objętości płaszczyzną odległą od środka bryły o fæl 

i określoną wektorem normalnym R. Skorzystajmy z przybliżenia 

podanego w dodatku Pa rusen wtedy 

V(I] wyw, ) = tuylar fashit S(r4+aB+tC- Iza] ) (8.2) 
8 ZA -A 


wykorzystując dwie równości; 


J6)- Gla) OG), (B3) 
(FB) = Fla) a) - Fla) Olz) (84) 


dla z») z, i F (z) oTi dz’ /f(z) jest dowolną funkcją całko- 
o 
walna/. Otrzyńujemy dla A,B,CfO /patrz dodatek A/ 


wzór: 


V(lzwlexcn)= (Ex AGABC| é 3 ar (zj > lz)” OLzje- fel) 
et (8. 5) 


gdzie Zijk" 1) JE €1) I B + (-1) ko „ Ponieważ V jest syme- 
tryczne w |Al, [B], [c], załóżmy że|A| 2 [B| > C 1 rozpiszmy (B.5) 
na poszczególne przypadki: 


y Iza| ZIAIŁIB|+|Q , (8.6) 
V=0; 
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2/ /trójkąt/ |AI+|B|-IC| 4 22] < [A|+]BI+IC| , 
Vee 
3/ gagy 14]- IBI + iq < [z2] Ś [A1+IB|- IC), 
Ve CAL IB] - [2421); 
4/ /pieciokat/ | 14)- I8|-Ic| |< 9 Ś [A|-1B] +le] 
V; 2 

V= jejapo | (A+ IB] + [C] -1z421)— 

— (IA +1B]-1c|- Jz) - 

- (IA-IBI+IC|- 12421)" Í 


ss [za]< |AI-IB|-Icl| 
a. /równoległobok/ [Ą| > IBl+[C] 


V= Vo /21A] 
b. /sześciokąt/ |A| < [Bl+lc] 
V= rega | (IA! + Bl+ le] 120) 
— (HAL + IB] IE] - Izl)? 
- (IA|-IBI+IC| - alf- 
- CIA] + [B) Hd = Zel Å | 


(8.3) 


(8.9) 


(B.3) 


(8.10) 


(B.n) 


Y nawiasach podaliśmy nazwy figur płaskich,będęcych przekrojami 


bryły, Jeśli A,BĘO a C=0 tylko przypadki 1,3 i 5a pozostają. 


Jeśli Ał0 a B=C=0,pozostają tylko przypadki 1 i 5a. 
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Dodatek C 


Współczynniki rozwinięcia 7 å af) „ ze-wzoru (6.33) /przypadek 


idealnego uporządkowania orientacyjnego/ 


2 ke) 2 , 7 
g = ©) - żę: ze ee (gp) f {A-Aa} 
ag? W(k*) = Aka (n*) y? wk) WK a 


— aż? Ay (8) BR BOR] / {1+2 Atl) Ma 


tt AW, (qe) ne), EN 
042 2'F z 2) JF 
aż -($) 20 €=0 39 R le-o = 
23 < 
ng le-o 


3 Å l p- (g2)* E AV (5) RARA - 
- 4 Abe i R) (o) - Aby (GP 
« WOKE) NAKE) = $ Aka (E) 4 Ree) (6) - 
- Ż My (1) 3 RUKE) (Q-A r GN? - 
+ KY CKE) GP — Ad) gt? AS / 
/ 12044 (ARK) +2 Abe (FINK) + 


+ Ary (7) mA) + 2 AV, (ER 4x9) | (2) 
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